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ABSTRACT 

The  spectral  analysis  technique  is  used  in  order  to  predict  the 
wave  induced  rolling  motions  of  a  ship  in  a  seaway.   Non-linear 
equations  of  motion  for  the  ship,  the  anti-rolling  tank,  and 
for  wave  excitations  are  developed.   Solution  techniques  for 
the  equations  are  examined,  and  a  linearized  form  is  adopted. 
A  design  method  using  a  digital  computer  program  based  upon 
numerical  solution  of  the  linearized  equations  is  proposed  and 
examined  for  a  number  of  different  representative  ships.   The 
design  method  is  demonstrated  to  be  capable  for  use  to  design 
passive  anti-rolling  tanks  for  ships  and  to  determine  feasi- 
bility of  activated  anti-rolling  tanks. 
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I.  Introduction 

1.  General  Discussion 

Successful  roll  stabilization  of  ships  by  means  of 
fluid  filled  Anti-Rolling  tanks  requires  a  thorough  analysis 
of  the  forces  and  moments  involved.   The  author  undertook 
this  task  with  a  goal  of  being  able  to  use  the  results  as  a 
design  tool  of  the  naval  architect.   As  such,  the  Anti-Rolling 
tank  design  method  derived  in  this  paper  is  intended  to  be 
a  single  step  in  the  larger  process  of  designing  the  whole 
ship. 

Modeling  the  dynamics  of  a  ship  with  an  Anti-Rolling 
tank  has  been  investigated  by  many  authors  of  whom  only  a 
small  percentage  are  represented  in  the  reference  list.   Sev- 
eral problems  have  prevented  their  analyses  from  being  complete 
or  from  possessing  sufficient  accuracy  in  the  range  of  interest, 

One  central  problem  is  that  rolling  is  a  resonance 
phenomenon.   The  ship  is  an  underdamped  dynamic  system  in  roll, 
one  which  magnifies  wave  excitations  to  a  certain  degree.   This 
may  be  proved  by  considering  that  it  is  well  established  for 
waves  in  the  open  ocean  to  have  a  maximum  waveslope  of  about 
16  degrees.   Waves  much  steeper  than  16  degrees  are  unstable 
and  break  up  by  forming  whitecaps  or  similar  turbulence. 
However,  rolls  of  2  0  or  3  0  degress  or  more  are  common  even  in 
moderate  seas.   Thus,  for  the  normal  range  of  ship  parameters 
and  dimensions,  a  resonant  condition  occurs  within  the  fre- 
quency range  of  wave  exciting  forces. 

To  better  understand  the  implications  of  being  in  a 


resonant  condition,  consider  the  ship  as  a  simple  second 
order  system  in  roll/  with  the  coordinate  <j>  representing 
roll  angle. 

(Kp-Ix)  "(J>  +  Kp$  +  K<t>4>  =  -  (Wave  excitations) 


The  natural  resonant  frequency  of  the  system  is  to  =^  /  Kcp 

(Kp-Ix) 


is  a)  =  / 


and  it  may  be  recognized  that  the  amplitude  of  the  roll  motion 
at  the  resonant  frequency  is  proportional  to    vKo'  (Kp-Ix)  /Kp 
Hence  the  magnitude  of  the  rolling  motion  at  resonance  will  be 
most  strongly  dependent  upon  the  damping  parameter  Kp.   Some 
analyses  of  ship  motions  have  assumed  a  linear  damping  para- 
meter for  the  ship,  i.e.  Kp  =  a  constant.   This  is  a  good 
approximation  for  small  motions  or  for  specific  conditions, 
but  in  the  general  case  we  expect  Kp  to  be  non-linear.   Since 
roll  damping  involves  viscous  friction,  we  expect  that  there 
will  be  a  Reynolds  number  dependency.   Also,  because  rolling 
creates  waves,  there  is  added  damping  due  to  the  radiation  of 
wave  disturbances  from  the  ship.   This  radiation  of  wave  energy 
has  been  shown  to  be  frequency  dependent  in  a  complicated 
fashion. 

Another  problem  is  that  rolling  motion  is  not  indepen- 
dent of  other  ship  motions  and  therefore  cannot  be  reasonably 
expressed  as  a  simple  second  order  system.   For  example,  a 
ship  at  rest  has  been  caused  to  assume  an  angle  of  roll  which 
is  different  from  its  usual  upright  position.   This  roll  angle 
causes  the  ship  to  have  an  unsymmetrical  underwater  shape, 
which  induces  a  righting  moment  (assuming  a  stable  ship) . 
Obviously,  any  pitching  motion  would  change  the  underwater  shape 


still  further,  unless  the  ship  were  prefectly  spherical, 
and  would  induce  a  change  in  the  value  of  the  righting  moment. 
It  can  be  shown  by  similar  arguments  that  roll  is  coupled  to 
all  of  the  other  ship  motions  and  displacements.   However, 
extensive  studies  have  shown  that  the  degree  of  coupling  is 
negligible  in  all  instances,  except  for  sway  rate  and  sway 
acceleration,  assuming  a  "normal"  surface  ship. 

This  paper  will  take  an  engineering  approach  towards 
analyzing  rolling  motion.   The  best  available  information  on 
how  to  model  non-linearities  will  be  used  in  a  workable 
system  which  allows  calculations  to  be  done  in  a  reasonable 
amount  of  time.   This  approach  was  taken  because  of  the 
iterative  nature  of  the  design  process  for  Anti-Rolling  tanks., 
which  in  turn  is  one  step  in  a  larger  iterative  design  process 
for  the  entire  ship.   Therefore  the  analytical  model  of  the 
Anti-Rolling  tank  had  to  be  flexible  enough  to  be  "turned 
inside  out"  into  a  systhesis  model  by  which  a  large  range  of 
parameters  could  be  investigated  by  some  search  technique. 
This  factor  was  the  primary  motivation  behind  taking  an  engin- 
eering approach  to  a  sufficient  design  rather  than  a  purely 
analytical  approach  towards  a  perfect  dynamic  model.   It  is 
hoped  that  this  paper  possesses  a  good  enough  balance  between 
numerical  estimation  and  reality  to  shed  a  little  more  light 
on  a  very  complicated  subject. 
2 .  System  of  Units 

The  English  system  of  units  is  used  throughout  this 
paper  and  in  the  computer  program.   The  following  English  units 


are  assumed  pertinent  to  all  applicable  quantities  unless 
other  units  are  indicated. 

Time:   seconds 

Distance:   feet 

Mass:   slugs 

Angle:   radians 

2 

Gravity  constant:  32.2  ft.  per  sec. 

Displacement:   long  tons  (2240  lbs) 
Force:   pounds 
Moment:   foot-pounds 


II.  Derivation  of  the  Equations  of  Motion 
1.  General  Discussion 

This  chapter  presents  a  detailed  derivation  of  the 
equations  of  motion  of  the  ship-tank  system.   The  form  of 
the  equations  is  highly  dependent  upon  the  coordinate  system 
chosen,  which  makes  it  important  to  understand  section  2  and 
figures  1  and  2  before  continuing  with  the  text.   For  those 
who  are  interested  only  in  the  results,  section  9  presents 
the  final  form  of  the  equations  and  section  10  gives  their 
non-dimensionalized  equivalent.   It  is  suggested  that  readers 
who  are  not  interested  in  the  derivations  should  still  read 
section  2  before  skipping  to  the  results.   This  will  ensure 
the  reader  has  a  reasonable  feeling  for  the  basic  nomenclature 
of  the  equations  before  facing  them. 

The  derivation  of  the  equations  of  motion  begins  with 
Newton's  second  law  of  motion,  action  equals  reaction.   All 
of  the  terms  which  make  up  the  reactive  forces  and  moments  are 
first  listed.   These  terms  are  each  broken  down  into  their 
dependent  variables.   An  expression  for  the  reactions  is  then 
formed  by  taking  the  first  partial  derivative  of  the  partic- 
ular force  or  moment  with  respect  to  each  of  it's  dependent 
variables.   If  the  indicated  partial  derivatives  were  all 
constants  the  reaction  terms  would  be  exactly  linear.   Each  of 
the  partial  derivatives  is  examined  for  this  property,  simpli- 
fications are  made  where  possible  and  in  several  cases  secondary 
parameters  are  introduced  to  further  define  particular  coeffi- 
cients.  Enough  documentation  of  relevant  assumptions  is  made 


to  be  able  to  change  particular  coefficient  values  whenever 
the  state  of  the  art  permits.   The  excitation  forces  and 
moments  are  handled  in  the  same  fashion  except  that  they 
have  independent  as  well  as  dependent  variables. 

The  last  section  of  this  chapter  gives  the  summarized 
results  as  a  set  of  simultaneous  equations  which  are  in  non- 
dimensional  form.   This  section  is  the  key  reference  in 
understanding  the  rest  of  this  paper.   This  is  because  the 
non-dimensional  coefficient  parameters  found  there  constitute 
a  kind  of  universal  set  which  not  only  appears  in  most  other 
literature  on  Anti-Rolling  tanks  (often  with  the  same  symbology) 
but  also  is  used  continuously  throughout  the  rest  of  the  paper. 
2.  Coordinate  System 

The  coordinate  axes  for  the  ship  are  through  the  center 

of  gravity  of  the  ship  and  lie  along  the  principal  axes  of 

inertia.   They  are  x,y,  and  z.   The  reference  coordinates 

with  respect  to  the  earth  are  x  ,  y  ,  and  z  ;  positive  z 

o   J  o       o  o 

direction  is  in  the  direction  of  gravity. 


Figure  1.   Ship-Oriented  Coordinates 

The  tank  coordinate,  n ,  is  relative  to  the  ship  and  is 
positive  in  the  direction  shown  on  Figure  2.   All  of  the  symbols 


shown  on  Figure  2  are  in  their  positive  sense. 

The  intent  of  this  figure  is  to  show  the  relative 
positions  of  points  G,  B,  and  0.   G  denotes  the  center  of 
gravity  of  the  ship  with  a  filled  Anti-Rolling  tank  (cross 
duct  blocked)  and  serves  as  the  reference  for  the  ship 
oriented  x-y-z  axes.   B,  the  center  of  buoyancy  according  to 
Archimede's  principle,  is  assumed  to  be  the  point  of  appli- 
cation of  the  wave  excitation  forces  and  moments.   0  is  the 
origin  of  the  tank  coordinate  n .   Also  shown  is  the  waveslope 
angle  \p   and  some  of  the  basic  Anti-Rolling  tank  geometry. 
Figure  A-l  in  appendix  A  gives  a  more  complete  picture  of  the 
Anti-Rolling  tank  by  itself. 
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3.   Reaction  Forces  and  Moments. 

For  the  ship-oriented  coordinates  a  force  and  moment 
balance  is  done.   K  denotes  a  roll  moment  about  the  x-axis, 
Y  denotes  the  force  in  the  y  direction,  and  E  denotes  the 
tank  moment  about  an  axis  parallel  to  the  x-axis. 

I    (Moments  in  <j>  direction)  =  K  =  I<j> 

E  (Forces  in  the  y  direction)  =  Y  =  mv 

For  the  tank  coordinate,  reference  (4)  uses  a  Lagrangian 
description  of  the  potential,  kinetic,  and  dissipation 
energies  associated  with  the  tank  to  arrive  at  its  equation 
of  motion.   For  the  sake  of  brevity,  only  the  results, 
which  are  well  known,  will  be  used  in  this  paper. 

Assuming  a  full  tank  with  blocked  ducts,  there  will 
be  no  interaction  between  the  tank  and  roll  or  sway.   In 
all  other  cases,  there  will  be  coupling  between  the  three 
coordinates.   This  leads  to  an  expression  for  K  and  Y  as 
follows: 

K  =  K(tank)  +  K (hydrodynamic)  +  K (excitations) 
Y  =  Y(tank)  +  Y (hydrodynamic)  +  Y (excitations) 

The  hydrodynamic  forces  and  moments  are  those  imparted 
to  the  body  of  the  ship  by  virtue  of  the  ship  motions 
occuring  within  and  at  the  surface  of  a  dense  fluid  (sea- 
water)  in  the  presence  of  gravity.   The  tank  forces  and 
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moments  are  caused  by  the  tank  fluid  shifting  its   weight 
transversely  in  the  ship.   Free  surface  waves  at  the  surface 
of  the  ocean  are  the  cause  of  excitation  forces  and  moments. 

This  separation  of  terms  in  the  force  and  moment 
equations  may  be  carried  further  by  noting  which  variables 
correspond  to  tank,  hydrodynamic  or  excitation  portions. 

K  or  Y(tank)  =  K  or  Y  (n,  n,  n) 

K  or  Y  (hydrodynamic)  =  K  or  Y  (<j>,  <p ,  § ,    v,  v) 

K  or  Y  (excitations)  =  K  or  Y  (ij>) 

The  usual  process  of  obtaining  a  first  order  model 
of  a  system  is  to  take  the  first  partial  derivatives  of  each 
quantity  (K  or  Y  in  this  case)  with  respect  to  each  of 
its  arguments.   After  following  the  above  procedure, 
combining  terms  where  possible,  and  separating  all  the  reactive 
forces  and  moments,  including  those  of  the  tank,  from  the 
excitations,  the  roll   and  sway  equations  are  obtained. 

(roll)   (KvD  +  Kv)v  +  [  (Kp  -  Ix)D2  +  KpD  +  K<j)](J>  +  K(tank)  = 

-K (excitations) 
(sway)   [  (Yv  -  m)D  +  Yv]v  +  YpD2  +  YpD  +  Y<j>)<f>  +  Y(tank)  = 
-Y (excitations) 
where:   D  =  d/dt 

Kv  =  3/(3v)K       Yv  =  3/(3v)Y 
Kv  -  3/ (3v)K       etc. 
Kp  =  3/(3p)K 
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Kp  =  3/(3p)K 
K<f>  =  3/(3<(»)K 

The  first  order  equation  of  .motion  for  the  Anti-Rolling 
tank,  as  found  in  reference  (4)  is: 

(tank)    pAQR2l2Dv  -  (S"D2  +  2g)  <J>  -  (S'D2  +  BTD/2pA0R2  + 
2g)nl  =  -EQ(t) 

S'  and  S"  are  geometry  factors  which  are  calculated 
in  Appendix  A. 

The  coupling  terms  between  roll,  sway,  and  the  tank  are: 

K(tank)  =  Knn  +  Knn  +  Knn 
Y(tank)  =  Ynn  +  Ynn  +  Ynn 

Since  the  coupling  terms  must  be  symmetric  in  any  energy 
conservative  system,  the  tank  equation  gives  us  directly 
the  indicated  partial  derivatives  in  the  roll  and  sway 
equations. 

Kn  =  -pAqR^1'  Yn  =  +2pAQR2 

Kn  =  0  Yn  =  0 

Kn  =  -2pgAQR2  Yn  =  0 

Using  these  results  and  taking  the  case  of  a  purely 
passive  Anti-Rolling  tank  (E  (t)  =  0)  the  linear  equations 
of  motion  for  the  ship-tank  system  are  formed.   The  result  is 
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a  fifth  order  differential  equation  which  is  linear  if  the 
assumption  is  made  that  the  coefficients  of  all  the  terms 
are  constant. 

(roll)    (KvD  +  Kv)v  +  [Kp  -  Ix)D2  +  KpD  +  K<j>]<J>  -  pAQR2  (S"D2  + 

2g)n  =  -KQ(t) 
(sway)   [  (Yv  -  m)D2  +  Yv]  v  +  (YpD2  +  YpD)  <f>  +  2pAQR2D2n  = 

-YQ(t) 
(tank)   2pAQR2Dv  -  pAQR2  (S"D2  +  2g)  <f>  -  pA0R2(S'D2  +  B^/ 

pAQR2D  +  2g)n  =  -EQ(t)  =  0 

In  order  to  add  a  control  system  to  the  tank  (mathematically) 
it  is  necessary  to  set  the  right  hand  side  of  the  tank 
equation,  -Eq (t) ,  equal  to  -A  Rp (t) ;   p(t)  is  the  pump 
pressure  developed  across  a  control  effector  in  the  tank 
duct,  with  a  (+)p  direction  in  the  (-)n  direction.   A 
separate  equation  for  the  dynamics  of  the  pump  must  be 
included  in  the  equations  of  motion.   A  simplified  form  of 
p(t)  which  should  be  sufficient  in  the  general  case  is  to 
consider   that  any  control  effector  can  be  characterized  in 
a  gross  sense  by  an  equivalent  first-order  time  lag.   This 
leads  to  the  equation: 

Tpp(t)  +  p(t)  =  pc(t) 

pc  (t)  =  f  (<j>,$,v,n,n,tM 

Pc(t)    =    fg1(j>   +  g24>   +   g3n   +  g4n   +   g5v  +  g^,     |pc|<   Pmax 

i-+Pmax'     1^1*    +-"+g6^l     >    Pmax 
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The  expression  for  pc (t)  was  chosen  based  on  linear 
regulator  theory  which  allows  feedback  of  all  the  state 
variables,  plus  the  excitation  ifr.   In  practice,  i>    is 
random  and  unmeasureable  (economically) ,  and  of  the  state 
variables,  the  roll  and  tank  variables  are  the  easiest  to 
measure  and  the  most  strongly  coupled.   This  reduces  the 
controller  equation  down  to: 

Pc(t)  =  gx*  +  g2$  +  g3n  +  g4n,  | pc |  <  pmax 
Pc(t)  =  ±Phiax'  1^1+  +  $2*  +  93^  +  94^1  1  Pmax 

A  useful  equation  which  must  be  considered  along  with 
the  above  set  of  state  equations  is  an  estimation  of  the 
power  requirements  of  the  control  effector.   Using  the 
relationship  that  power  is  the  product  of  pressure  differ- 
ential and  volume  flow,  and  taking  into  account  pump  efficiency, 
power  requirements  are  calculated  by: 


Pump  Power  =  (1/e  )p(t)A0Rn 

ir 

where:   e^  -  0.5  -  0.8 


4.   Simplifications  of  the  Hydrodynamic  Coefficients. 

Assume  that  the  roll  and  sway  motions  of  the  ship 
as  seen  at  the  center  of  buoyance  are  uncoupled;  that  is, 
a  roll  motion  about  the  center  of  buoyancy  induces  no 
swaying  and  likewise  a  sway  motion  at  the  center  of  buoyancy 
induces  no  rolling  there.   It  may  be  easily  imagined  that 
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this  is  nearly  true.   Using  this  basic  assumption,  many 
simplifications  in  the  terms  of  the  hydrodynamic  coefficients 
are  possible.   For  example ,  a  roll  motion  about  the  x-axis, 
which  passes  through  the  center  of  gravity,  will  be  seen  at 
the  center  of  buoyancy  as  a  roll  motion  plus  a  sway  motion. 
The  hydrodynamic  resistance  to  sway  at  the  center  of  buoyancy 
is  felt  as  a  sway  force,  but  does  not  couple  back  into  roll, 
with  this  assumption.    Using  the  (')  notation  to  indicate 
that  a  particular  coefficient  is  taken  at  the  center  of 
buoyancy  instead  of  being  referenced  to  the  center  of  gravity 
(which  is  the  center  of  the  roll  and  sway  axes)  the  values 
of  the  following  hydrodynamic  coefficients  are  determined. 

Y'v  =  Yv  Yp=Y'p=  -ZftYv 

Y'v  =  Yv  Yp  =  Y'p  =  -zhYv 

K'<|>  =  K<|>  Kv  =  K'v  =  -zhYv 

Kv  =  K'v  =  -z,Yv 
h 

Introduction  of  the  above  relationships  shows  which 
hydrodynamic  coefficients  must  be  known  as  a  minimum  in 
order  to  solve  the  equations  of  motion  of  the  ship-tank 
system.   Notice  that  these  coefficients  are  all  defined 
as  being  determined  with  respect  to  the  x-y-z  coordinate 
system  about  the  center  of  gravity  of  the  ship. 

K<(>  =  Static  righting  moment  characteristic  of  the  ship. 
Kp  =  Linear  damping  in  roll. 
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Kp  =  Added  inertia  of  the  ship  in  roll  to  accomodate  water 
flow. 

Yv  =  Linear  damping  in  sway. 

Yv  =  Added  inertia  of  the  ship  in  sway  to  accomodate 
water  flow. 


5.   Tank  Coefficients. 

The  parameters  which  determine  the  coefficients  of  the 
tank  variables  in  the  linear  equations  of  motion  are: 

pAqr2  =  Characteristic  of  tank  size. 

S"  =  Function  of  tank  geometry  =   /  (d/R)dl 

S'  =  Function  of  tank  geometry  =   /  (Ar)/a)dl 

Bt  =  Tank  linear  damping. 

S'  and  S"  are  the  weighted  moment  arm  and  sectional  area, 

respectively,  of  the  tank  fluid.   They  are  familiar  terms 

2 
in  literature  on  hydraulic  engineering.   PA0R   has  the  units 

of  mass  times  a  moment  arm  and  is  a  characteristic  tank 

size  or  inertia  term. 

The  above  tank  parameters  are  all  free  to  be  chosen  . 

However,  it  can  be  noticed  that  S"  and  S1  fall  within  the 

following  range  of  possible  values. 

1  1   s'  1  °° 
-1  <  S"  <    1 

lBtl  1  O.OSCv^gS1"  pA0R2) 
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For  the  case  of  S'  getting  very  large,  the  tank  dynamics  slow 
down  due  to  the  increased  fluid  inertia.    The  extreme  case 
is  the  equivalent  of  the  tank  being  just  an  additional 
weight  aboard  ship.   The  sign  change  of  S"  is  accomplished 
by  changing  the  center  of  gravity  of  the  tank  relative 
vertically  to  the  center  of  gravity  of  the  ship.   A  negative 
sign  indicates  a  high  tank,  a  zero  value  occurs  near  the 
center  of  gravity  of  the  ship,  and  a  positive  sign  is  caused 
by  a  low  tank. 

The  damping  parameter  of  the  tank,  B, ,  is  free  to  be 
chosen  as  long  as  it  exceeds  the  approximate  minimum  shown 
above,  say  about  5  percent  of  critical  damping.   Later 
analyses  will  show  that  a  good  passive  tank  will  incorporate 
far  more  damping  than  this.   Otherwise,  undesirable  secondary 
resonance  peaks  appear  in  the  frequency  response  of  the 
ship-tank  system.   Therefore,  for  practical  purposes  B,  may 
be  considered  a  free  parameter  with  the  stipulation  that 
in  the  final  design  of  the  tank  suitable  damping  devices  such 
as  sharp  corners,  baffles,  etc.  are  used  to  achieve  the 
desired  degree  of  energy  dissipation.   B,  may  be  calculated 
or  predicted  accurately  enough  with  standard  hydraulic 
engineering  'pipe  friction'  tables. 

6.   Effector  Pump  Parameters. 

The  maximum  efficient  effector  size  can  be  easily 
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estimated  for  any  particular  tank  geometry,  since  the  effector 
pump  is  installed  in  the  cross  ducting.   This  will  enable 
the  choice  of  a  particular  sized  pump  which  we  have  chosen 
to  characterize  by  its  maximum  pressure,  Pmavr  and  lag  time, 
T  .   There  is  no  particular  requirement  for  Pmaxr  except 
that  pmax  =  (Power,  max. )  ep/A0Rri  max.   For  T   however  there 
is  the  requirement  that  the  pump  time  lag  must  be  less  than 
the  natural  resonant  frequency  of  the  tank  for  the  pump 
to  have  much  useful  effect.   A  pump  with  more  lag  will 
introduce  instabilities  into  the  control  system. 

7.   Non-Line ari ties. 

The  hydrodynamic  coefficients  obviously  all  depend  upon 
ship  speed  because  they  are  caused  by  changing  the  flow 
field  of  the  water  about  the  ship.   Of  these,  the  damping 
parameters,  Yv  and  Kp,  are  the  most  speed  dependent  because 
they  involve  the  dynamic  lift  in  the  horizontal  plane  generated 
by  the  ship  hull  serving  as  an  inefficient  lifting  body. 
An  angle  of  attack  is  imparted  to  the  mean  water  flow  by 
the  relative  motion  between  the  ship  and  the  water  in  the 
horizontal  plane  determined  by  the  vector  addition  of  sway 
rate  and  forward  speed.   The  resultant  lift  force  is  felt 
as  a  resisting  force  in  the  sway  direction.   The  static 
righting  moment  coefficient,  K(j>,  is  the  least  speed  dependent 
because  the  only  change  from  zero  speed  conditions  is  that  the 
constant  water  pressure  surfaces,  such  as  the  air-water  inter- 
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face,  are  no  longer  planar,  but  are  slightly  folded  into  a  wave 
pattern.   Based  on  these  arguments,  and  with  no  further  proof, 
it  was  decided  to  adopt  Kp  and  Yv  as  ship  speed  dependent,  and 
the  other  hydrodynamic  coefficients  as  speed  independent. 

Because  the  magnitude  of  the  roll  angle  at  resonant  fre- 
quenty  is  much  larger  than  its  average  value,  the  value  of  K^ 
does  not  remain  reasonably  constant,  say  within  5-10%,  within 
its  expected  range.   Since  the  well  known  curve  of  static  K 
versus  0  is  anti-symmetric,  it  was  decided  that  it  would  be  suf- 
ficient to  adopt  an  odd  polynomial  expression  for  K^ ,  as  was 
done  in  reference  (7).   The  Kg$0  term  is  replaced  by: 

The  K^3  coefficient  may  be  simply  evaluated  by  noting  that 
the  value  of  K  is  zero  when  0  is  at  the  maximum  static  equilibrium 
angle,  0r  .   This  fixes  the  value  of  K^3  as  -K0/0R2  .   This  can 
be  even  further  simplified  by  approximately  0R  by  about  80°  or 
1.4  radians.   This  enables  the  use  of  K^3  =  -K^/2  if  better  in- 
formation is  lacking. 

Part  of  the.  energy  dissipation  caused  by  roll  motion  is 
due  to  viscous  friction  resisting  the  relative  motion  of  the 
skin  of  the  ship  to  the  water.   Since  viscous  friction  is  Rey- 
nolds Number  dependent,  it  must  be  roll  velocity  dependent.   This 
component  of  roll  damping  can  be  calculated  by  standard  hydraulic 
engineering  'friction  factor'  tables  to  include  skin  friction 
and  the  eddy-making  resistance  of  the  bilge  keels.   However, 
this  usage  precludes  us  from  making  the  odd  polynomial 
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approximation  because  the  standard  tables  are  all  of  the  form 

o 
l/2p (area) (velocity  ) (friction  factor).   This  forces  the 

adoption  of  replacing  Kpp  by: 

KpP  +  KjppPIp 
The  absolute  value  sign  is  necessary  to  make  the  symmetric  p^ 
function  an  anti-symmetric  function.   Since  the  KrJpjPiP1  term  is 
primarily  due  to  transverse  flow  meeting  the  tremendous  flow 
obstruction  of  the  bilge  keels,  and  not  for  most  cases  much 
influenced  by  forward  speed,  it  was  decided  to  let  this  quantity 
remain  independent  of  ship  speed. 

Another  form  of  energy  dissipation  by  rolling  or  swaying 
motion  is  found  in  the  formation  of  surface  waves  which  radiate 
outward  from  the  ship.   This  energy  dissipation  has  been  esti- 
mated, using  potential  flow  theory,  for  simple  ellipsoidal 
bodies  which  resemble  the  underwater  form  of  ships,  and  has 
been  found  to  depend  on  ship  speed  and  the  frequency  of  motion 
(assuming  harmonic  motion) .   As  would  be  expected,  sway  gener- 
ates larger  waves  for  an  average  hull  form  and  therefore  this 
effect  is  more  pronounced  for  sway.   Thus,  it  was  decided  to 
approximate  sway  damping  as: 

(Yv1(uj,U)  +  Yv2(U))v; 
where  Yv,  is  the  component  associated  with  wave  formation  and 
YV2  is  the  component  associated  with  the  hull  acting  as  a  lift- 
ing surface  in  the  horizontal  plane.   Roll  damping  may  be  well 
enough  approximated  by: 
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KpP  +  Kp^ptp); 

where  Kp  is  calculable  by  strip  theory,  or  may  be  estimated 
from  Appendix  B  as  about  .05  times  critical  damping.   Kpjpi  can 
be  estimated  from  friction  factor  tables  at  an  average  Reynolds 
number.   If  the  above  seems  inconsistent,  it  should  be  recalled 
that  roll  damping  is  much  less  than  the  amount  of  sway  damping 
that  is  coupled  into  roll  by  the  equations  of  motion.   Thus, 
the  extra  care  is  taken  in  the  estimation  of  sway  damping  vice 
roll  damping. 

To  determine  the  above  coefficients  exactly  requires  that 
extensive  model  testing  be  conducted.   Since  it  is  envisioned 
that  the  method  proposed  by  this  paper  is  done  at  a  stage  of 
design  prior  to  any  such  detailed  knowledge  of  the  ship,  the 
coefficients  must  be  arrived  at  either  by  empirical  data 
derived  from  previous  similar  designs,  or  else  by  estimation 
from  their  theoretical  values  for  similar  geometric  forms.   An 
appendix  to  this  paper  gives  detailed  examples  of  how  to  make 
these  estimations. 

.As  far  as  the  tank  non-linearities  are  concerned,  the 
restriction  on  the  maximum  tank  angle,  n ,  by  the  boundaries  of 
the  tank  is  the  most  important  of  the  non-linearities.   The 
coupling  parameters  Sand  S"  remain  almost  constant  for  any  tank 
geometry  as  long  as  centerline  symmetry  is  preserved.   However, 
a  saturation  effect  is  observed  whenever  the  fluid  in  the  tank 
reaches  the  tank  top  because  at  that  point  no  more  fluid  transfer 
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athwartships  is  allowed  so  that  the  tank  moment  reaches  a  maxi- 
mum.  To  account  for  saturation  in  the  equations  of  motion  there 
are  two  methods.   An  equivalent  linear  'spring  constant1  may  be 
used;  to  do  this  change  S'  such  that  the  same  energy  storage 
for  a  given  tank  angle  is  achieved  for  a  fictitious  unsaturated 
tank  as  is  realized  in  the  real  tank.   Another  method  is  to 

change  the  equations  of  motion  at  the  point  that  the  tank  angle 

•    •• 

saturates.   After  saturation,  consider  that  n  =  n  =  0 .   The 

normal  tank  dynamics  resume  again  when  the  fluid  can  run  'down- 
hill'.  This  condition  is  met  in  the  static  case  for  n  +  $    >    0. 
In  the  dynamic  case,  the  dynamic  head  must  be  included  along 
with  the  pump  pressure,  if  present.   The  condition  for  resumption 
of  normal  tank  dynamics  then  is: 

2PAqr2dv  -  pA0R2(S"D2  +  2g)0  -  2pA0R2gn 
+  A0Rp(t)  >  0. 

The  tank  damping  parameter,  B^,  of  course  has  non-lineari- 
ties associated  with  it  due  to  the  nature  of  the  damping  devices 
used  in  the  tank.   However,  as  mentioned  earlier,  because  we 
are  interested  in  tank  design,  B-f-  can  be  considered  a  linear 
coefficient  independent  of  the  other  variables  with  the  stipu- 
lation that  the  amount  of  linear  damping  must  be  later  designed 
into  the  tank  by  some  means.   It  is  reasonable  to  expect  that  an 
amount  of  damping  in  the  same  order  of  magnitude  as  critical 
damping  is  easily  achievable. 
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8 .   Excitations 

In  order  to  keep  the  exciting  function  one-dimensional, 
it  was  arbitrarily  decided  to  reference  all  the  excitations  to 
the  waveslope.   This  was  mostly  done  to  be  able  to  compare 
results  with  previous  work.   The  waveslope  used  was  not  the 
surface  waveslope,  but  an  average  value  which  the  ship  senses 
(remembering  that  surface  waves  are  a  local  phenomenon  which 
decay  exponentially  with  depth) .   This  average  value  is  often 
taken  at  the  half-draft  of  the  ship,  but  to  be  consistent  with 
hydrostatics  as  well  as  hydrodynamics,  I  chose  to  use  an  aver- 
age value  as  the  value  at  the  center  of  hydrostatic  buoyancy, 
which  is  nearly  the  half-draft.   Also,  this  places  the  ex- 
citations at  the  center  of  buoyancy,  and  in  accordance  with  a 
previous  assumption  about  the  non-coupling  of  sway  and  roll  at 
this  location,  this  neatly  separates  roll  and  sway  excitations 

It  is  now  necessary  to  review  linear  wave  theory  for 
simple  harmonic  ocean  waves  in  deep  water.   A  later  assumption 
to  be  made  is  that  the  real  ocean  is  built  from  a  spectrum  of 
these  simple  harmonic  waves. 
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Figure  3.   Wave  Profile 
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Figure  3  is  arrived  at  from  the  equation  of  surface  ele- 
vation of  an  harmonic  water  wave,  given  in  reference  3,  which 
corresponds  to  a  surface  of  constant  (atmospheric)  pressure. 
Horizontal  wave  velocity,  acceleration,  and  waveslope  are  found 
by  taking  the  appropriate  partial  derivatives. 
e  =  ea  cos (ks-wt) 
vw  =  ae/9t  +  90°  !  =-toea  cos(ks-wt) 
vw  =  3vw/5t  =-w2ea  sin(ks-cot) 
\l>   =    de/ds   =  -kea  sin(ks-tot) 

An  important  factor  to  notice  from  the  above  math  is  that 
at  any  given  point  in  the  wave  \\>   and  vw  are  in  phase  with  each 
other  but  are  out  of  phase  with  the  wave  horizontal  velocity  vw 

For  simplicity,  the  wave  variables  may  be  expressed  in 
complex  form,  with  the  implication  that  the  quantity  indicated 
is  the  real  part  of  the  complex  expression. 

e  =  ea  exp  (i  (ks-tot)  ) 

•  2 

vw  =  -we;   vw  =  id)  e;  ip  =  ike 

vw/^  ~  ~w/ik  =  ig/w 
vw/>  =  u2/k  =  g 

The  last  two  relationships  fix  vw  and  vw  as  functions  of  i> 
Assuming  that  the  excitations  occur  at  the  center  of 
buoyancy,  and  using  the  same  prescription  as  before: 
KQ(t)  =  K'^  +  K'vwvw  +  K'^wvw 

Y0(t)  =  Y'^  +  Y'vwvw  +  Y';wvw 
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There  is  an  unrealistic  assumption  in  the  above  derivation, 
however,  in  that  the  horizontal  wave  motions  are  not  impinging 
on  a  static  ship.   They  are  actually  influencing  a  moving  ship, 
which  suggests  that  for  the  equation  to  be  valid,  vw  and  vw 
must  be  replaced  by  vw  -  v1  and  vw  -  v,  respectively. 

By  making  the  assumption  that  the  hydrodynamic  coefficients 
of  the  water  moving  past  the  ship  are  the  same  as  for  the  ship 
moving  past  the  water,  and  by  adopting  the  classical  assumption 
of  Froudethat  Kl  =  K^,  all  of  the  coefficients  are  determined. 
Notice  that  Yl  vanishes  due  to  a  previous  assumption,  but  that 
K'vw  and  K1^    cannot  vanish  unless  the  moment  arm,  z^  ,  between 
the  center  of  buoyancy  and  the  center  of  gravity  (the  reference 
for  the  x-y-z  axes)  vanishes.   With  all  of  the  above  assumptions 
and  conditions  included,  the  excitations  may  now  be  found. 

V  =  v  -  zhP 
v1  =  v  -  z^p 
K0(t)  =  Krfip   +  zhYv(vw-v  +  zhp)  +  zhYv(vw-v  +  zhp) 
Yo(t)  =  Yv(vw-v  +  zhp)  -  Yv  (vw  -  v  +  zhp) 

All  of  the  terms  in  the  excitations,  except  for  those  in- 
volving i|),vw(i|/)  and  vw(i|>)/  may  be  moved  to  the  left-hand  side 
of  the  equations  of  motion  when  introducing  this  expression  for 
the  excitations  into  the  equations  of  motion.   When  combined 
with  like  terms,  this  doubles  the  value  of  some  of  the  hydro- 
dynamic  coefficients,  and  introduces  coupling  terms  into  some 
others.   Also,  the  excitation  expression  added  two  new  parameters 
instead  of  the  one  desired — the  waveslope,  and  the  phasing 
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difference  between  waveslope  and  horizontal  wave  motions. 

9 .   Equations  of  Motion 

The  results  of  the  previous  sections  are  summarized  to 
completely  form  the  equations  of  motion  for  an  anti-rolling 
tank. 

(roll)       [(Kp-Ix)D2  +  KpD  +  K^]0    -  zh[YvD  +  Yv1{u,\J)     +   Yv2(U)]v 
-  pA0R2(S"D2  +  2g)n  +  Kpyplp;  -  K^3  =  -KQ(t) 

(sway)      zh[YvD2  +  (Yv1(to/U)  +  Yv2(U))D]jz5  +  [(Yv-m)D 
+  Yvx(w,U)  +  Yv2(U)]v  +  2pAQR2D2n  =  -Y0(t) 

(tank)      -pA0R2(S"D2  +  2g)^  +  2pAQR2Dv 

-pA0R2(S,D2  +  (Bt/pAoR2)D  +  2g) n  =  0;  or-  A0Rp(t) 

(effector)  (TpD  +  l)p(t)  =    (q±    +   g2D)^  +  (g3  +  g4D)v, 

|ZgiXi|<  Pmax  +  Pmax,  otherwise 

(excitations)   KG(t)  =  -Kfity   +   zh (Yv^  +  Yv2) (vw  -  v  +  zhp) 

•    •  •  • 

+  zhYv(vw  -  v  +  zhp) 
Y0(t)  =  ~(Yvx  +  Yv2) (vw  -  v  +  zhp)  -  Yv(vw  -  v  + 

zhp) 

•       •  • 

where:    n  =  n  =    0and|n|  =  nmax,  whenever  |n|exceeds  nmax.; 
the  condition  for  returning  to  n  and  n  free  is 
whenever: 
2pAoR  Dv  _  PAoR2[S..D2  +  2q}0    -    2pgAoR2n  +  AoRp(t)  >  0; 
with  the  initial  condition  n  =  +  nmax  and  n  =  0  . 


26 


These  equations  may  be  linearized  by  dropping  the  terms 
involving  o7   and  pjpi  ,  and  the  restriction  on  tank  angle  and 
pump  pressure.   A  more  convenient  way  to  express  these  equations 
is  to  non-dimensionalize  them;  the  results  in  non-dimensional 
form  represent  a  family  of  solutions  which  help  to  give  an  in- 
sight of  the  general  design  problem. 

10 .   Non-Dimensional  Equations  of  Motion 

Because  further  analysis  in  this  paper  uses  mostly  the 
linearized  version  of  the  equations,  only  the  linear  terms  are 
given  below.   A  paragraph  at  the  end  of  this  section  gives  the 
non-linear  terms  as  well. 

Dimensionless  variables  were  chosen  as  follows: 
t  =  wst 
v*   =  v/wsB 
p*   =  p(t)2pgR 

p7  ,  n  ,  \\)    in  radians 
D   =  d(   )/dx 
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The  linearized  equations,  with  complex  notation  for  the 
excitations,  are: 


(roll) 


(sway) 


[(1    +    Ab2Y!)D2    +     (Ks    +    Ab2Ky)D    +    1] 0 
-    2Ab[Y]_D    +    Ky]v*    +    At(D2/yst2    +    1)  n 
=     [  (1    -    AbAy)     -    iAbA0]     ^>(t) 

-2Ab[YiD2    +    KyD]0    +     (y2D    +    2Ky)     v* 
-    AtD2/    pyt2n      =    (Ay   +    iAo)     ^  ( t ) 


(tank) 


At(D2/pst2    +    D0    -    AtD/yyt2    v 

t 


+    At(D2/yt2    +    KtD/yt    +    Dn    =    0;    or    At    p*  (t) 


(effector)         2pgR(wsTp   +    Dp*   =    (gi   +   wsg2D)^   +    (93   +   03S94D)  ri 


Coupling  constants 

At  =-2pgAoR2/K^ 

Ab  =  zh/B 

Ay  =  YvgB/K^ 

Ao    =     (Yvi    +    Yv2)gB/oaoK0 

Inertia  constants 

A!  =  B2Yv/(Kp  -  Ix) 

A2  =  B2(2Yv  -  m)/(Kp  -  Ix) 

Spring  constant  (ship) 
K^  =  -mgGM 


Damping  constants 


KS    =   Kp/~\/K^(Kp    -    lx) 
Ky    =    B2YvVK^(Kp    -    Ix) 
Kt    =    Bt/pAoR2  V29    s' 

Frequency   constants 
yt   =   ait/us 
Pst   =   ust/ws 
pyt    =    Vg7B/u)S 
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In  order  to  preserve  the  non-linear  equations,  the  follow- 
ing terms  need  be  added  to  the  indicated  equation,  and  also  the 
restriction  on  tank  angle  and  pump  pressure  applies. 


(roll)      Kqp;pl  +  K3^3 


where:   Kq  =  (coS2  Kp,pl)  /K^ 

K3  =  1/2  ,  or  1/0r2 

There  is  one  great  idiosyncracy  of  these  equations  which 
would  pass  unnoticed  were  it  not  for  the  presence  of  the  co0  term 
in  the  wave  velocity  -  sway  coupling  constant,  A0  .   Here,  oj0 
refers  to  the  absolute  frequency  of  the  waves  relative  to  the 
earth.   The  rest  of  the  dynamics  will  occur  at  a  steady-rate 
frequency  which  is  the  encounter  frequency  that  the  ship  senses 
in  the  waves.   This  encounter  frequency  is  a  function  of  both 
ship  speed  and  relative  heading  angle,  and  is  further  examined 
in  the  appendix  on  excitation  spectra.   The  peculiarity  is  that 
the  ship  is  sensitive  to  both  relative  and  absolute  excitation 
frequencies,  which  is  a  direct  result  of  the  choice  of  sway 
excitations . 
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III.   Excitation  Spectrum 

1 .   General  Discussion 

The  preceeding  chapter  completely  defined  the  anti- 
rolling  tank  problem  with  the  exception  of  the  single  inde- 
pendent variable  4>    (time)  .   This  chapter  is  dedicated  towards 
obtaining  a  realistic  expression  for  ip .   Ocean  waves  are 
essentially  a  random  process.   To  account  for  the  statistics 
of  the  wave  process  it  was  necessary  to  adopt  an  accepted 
empirical  spectral  (implying  frequency)  distribution  function 
which  had  the  disadvantage  of  being  an  elevation  distribution 
instead  of  a  slope  distribution.   The  main  part  of  the  text 
of  the  chapter  deals  with  making  corrections  and  additions 
to  this  distribution  function. 

It  is  first  changed  into  a  waveslope  distribution  by  use 
of  linear  gravity  wave  theory.   Correction  factors  are  derived 
for  the  amount  of  submergence  of  the  ship's  hull,  for  the 
relative  heading  of  the  ship  with  respect  to  wave  progression, 
for  the  wave  encounter  frequency  shift--akin  to  the  familiar 
doppler  shift--and  for  the  shape  of  the  ship's  hull.   Because 
of  the  highly  speculative  nature  of  some  of  the  assumptions, 
particular  attention  is  paid  to  their  expected  range  of 
validity.   All  of  the  factors  are  then  combined  to  form  a 
complicated  expression  for  the  excitations  which  can  be 
inserted  into  the  equations  of  motion  to  complete  the  model 
of  the  ship-tank  system  in  ocean  waves. 
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This  strategy  was  adopted  because  it  lumped  as  many  of 
the  parameters  which  are  independent  of  the  anti-rolling  tank 
design  as  possible  into  a  single  expression  for  excitations. 
This  could  then  be  used  as  an  input  forcing  function.   The 
resultant  form  of  ip  ■  (time)  had  the  peculiar  but  advantageous 
properties  of  being  "steady  state",  i.e.  a  combination  of 
regular  harmonic  waves,  and  linear  while  only  requiring  that 
the  various  auxiliary  component  parameters  be  analytic.   Thus 
the  final  derived  expression  for  the  excitation  spectrum  ip 
included  most  of  the  design  independent  parameters  as  well  as 
the  statistical  distribution  of  waveslope. 

2 .   Waveslope  Distribution 

In  order  to  describe  the  wave  excitations  of  the  sea 
surface  it  was  decided  to  use  the  empirically  derived  Pierson- 
Moskowitz  wave  amplitutde  spectrum.   This  is  a  standard 
analytical  formulation  which  seems  to  adequately  describe  the 
real  life  wave  amplitude  spectrum  for  fully  developed  seas. 
The  statistical  distribution  of  wave  amplitude  over  the  band 

of  wave  frequency  (0  to  °°)  is  assumed  to  depend  on  just  one 

.  .  1/3 

parameter,  significant  waveheight,   H    .   Significant  wave- 
height  is  commonly  used  as  a  measure  of  sea  state.   To  be 

2 
more  specific,  e  /6a),  wave  amplitude  squared  per  unit  fre- 
quency bandwidth,  is  a  function  of  frequency  only  for  any 
given  sea  state.   For  a  more  complete  description  of  this 
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spectral  density  function,  its  implications  and  proven  uses, 
consult  references  3  and  11. 


e2/6w  =  0.0081  g2to  5  exp(-6/w4) 


6  =  0.0324  g2/(H1/3)2 


2 
g  =  32.2  ft/sec. 


Using  linear  wave  theory,  it  was  previously  shown  that 


\\\)\    =  |9e/8s|  =   ke 


2 

where :     k  =  w  /g 


Making  the  appropriate  substitutions,  the  spectral  den- 
sity function  for  waveslope  is  found. 


i2/*     v2o   /  ^      -2-1  (-3/u)4) 

ty   /6oo  =  k  S„  (u)  =  ctg   <*>  e 


The  mean  square  value  of  the  waveslope  spectrum  is  found 
by  integrating  the  spectral  density  function  over  the  entire 
frequency  range. 


—2       oo   2 

ip    =   /   k  Srr(w)dw 
o      £5 
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The  root  mean  square  value  (RMS)  of  wave  slope  is  found 
by  \p  =jy2*      By  the  statistics  of  a  random  process  with  an 

assumed  Rayleigh  distribution,  the  average  of  the  1/N  highest 
values  of  ty   are  given  in  reference  1 : 


^1/N  =  f  (1/N)  •  i> 
r         '      rrms 


1/N  f  (1/N) 

1/1  1.25 

1/3  2.00 

1/10  2.55 

1/100  3.34 

1/106  5.27 


It  should  be  pointed  out  that  the  above  numbers  refer  to 
the  waveslope  angle  measured  from  the  horizontal,  and  not  to 
the  "crest  to  trough"  condition  usually  associated  with 
waveheights . 

Following  the  convention  that  "significant"  means  the 
average  of  the  1/3  highest  values  of  the  process,  the  signi- 
ficant waveheight  is: 

r         rrms 

This  represents  an  acceptable  peak  value,  for  engineer- 
ing purposes,  of  the  random  waveslope  process.   Notice  that 
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the  mean  square  value  has  a  different  relationship  in  this 
random  process  than  for  the  deterministic  case  in  which  peak 
value  is  /2  times  the  RMS  value. 

3.    Depth  Correction 

The  surface  elevation  of  the  sea  is  a  constant  pressure 
surface.   This  pressure  decays  exponentially  with  depth  until 
at  a  certain  depth,  the  surface  waves  are  not  felt.   By  linear 
theory,  this  exponential  decay  factor  is  e        /  where  z 
is  the  depth  from  mean  surface  elevation. 

Since  the  center  of  buoyancy  of  the  ship  is  the  point 
through  which  wave  excitations  are  induced  in  the  static  case, 
and  for  small  amplitude  excitations  (i.e.  nearly  static  case), 
a  reasonable  approach  towards  excitations  is  to  apply  them 
at  the  center  of  buoyancy.   The  only  alternative  to  this  is 
to  carry  out  the  integration  of  the  pressure  distribution 
over  the  entire  immersed  area  of  the  hull.   By  the  above 
reasoning,  the  spectral  density  function  which  the  ship  feels 
may  be  approximated  by: 


!-2tbuVg  s  {B)J 


where  t,  is  the  draft  to  the  center  of  buoyancy 


This  correction  factor  for  depth  is  shown  in  figure  4. 
Notice  that  for  super  ships,  with  drafts  of  80+  ft  and  t,  of 
about  40  ft,  waves  are  practically  not  felt  in  roll 
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excitation.   For  larger  ships,  say  20,000+  tons,  all  wave 
lengths  less  than  about  50-100  ft.  are  unfelt,  and  are  non- 
interactive  with  the  ship. 


20,000' 
Figure  4. 


200'       50.6* 
Depth  Correction 


22.5" 


CO 


X 


4 .   Encounter  Frequency 

Since  the  ship  forward  speed  and  wave  advance  are  both 
vector  quantities,  relative  encounter  frequency  of  the  ship 
with  the  waves  changes  with  orientation.   Denote  heading 
angle  of  the  ship  into  the  waves  by  3  (3  is  180°  in  head 
seas).   The  applicable  relationship  is: 


w   =  u)  (1  -  (u  U  cos  3/g) 
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5.  Transverse  Waveslope 

The  equations  of  motion  were  developed  for  roll  and  sway 
coordinates  which  have  entirely  a  transverse  sense  with 
respect  to  the  ship.   It  can  be  shown  that  for  waves  progres- 
sing at  some  other  direction  than  transversely  to  the  ship 
the  effective  waveslope  in  the  transverse  direction  is: 

ij>(3)  =  sin  (3);   3  is  90°  for  beam  seas. 

6.  Region  of  Validity  of  Excitation  Spectrum 

Combining  all  the  corrections  to  the  original  waveslope 
spectral  density  function,  an  expanded  function  is  formed. 


ip2/6w(u),  tb,  3)   =   sin2  (3)  e  2tbW  /g  S   (w) 


The  expanded  form  of  spectral  density  is  plotted  on 
Figure  6  for  beam  seas  (sin (3)  =  1)  for  the  normal  range  of 
expected  sea  states.   A  glance  at  the  figures  shows  that  for 
all  ships  with  a  draft  of  5  feet  or  greater,  the  desired 
range  of  integration  of  S^  is  0.2  <  u  <  3.0  to  establish 
ip    .   The  limits  on  oo  impose  the  corresponding  wavelength  to 
fall  within  20  5  A  <_  20,000  feet. 

To  determine  the  range  of  3  possible,  a  visualization  of 
the  ship  in  long  crested  waves  is  helpful. 
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Figure  5.   Ship  in  Long  Crested  Waves. 
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Significant  Wave  Height  =  30  ft 


oo 


Significant  Wave  Height 
Beam  Seas 


=  15  ft 


tb=0 
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7&u> 


Significant  Wave  Height 
Beam  Seas 


=  5  ft 


II 


1/3=5' 


tb=0 


3     <^ 


Figure  6.   Excitation  Spectral  Density  in  Long-Crested  Waves 
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It  seems  reasonable  to  assume  that  the  excitation  at  the 
center  of  buyoyancy  is  equivalent  to  that  for  the  whole  ship 
as  long  as  the  value  of  L|cos($) |  is  less  than  a  half- 
wavelength.   Another  consideration  is  that  the  roll  excita- 
tion, the  more  significant  of  the  excitations,  is  propor- 
tional to  Ka.   A  look  at  the  terms  in  Ka  is  enlightening. 


Ku  =  pg (Immersed  Volume)  GM 

K.  =   (a  constant)  GM 

BM  =   (2nd  moment  of  waterplane  area  about  x-axis) 

(Immersed  Volume) 

GM  =   BM  +  KB  -  KG 


The  above  implies  that  K   is  proportional  to  the  2nd 

moment  of  the  waterplane  area.   It  can  be  shown  that  for  a 

3 
rectangle,  the  2nd  moment  is  LB  /12.   Therefore,  it  can  be 

3 

seen  that  K   is  proportional  to  (beam)  . 

Since  the  beam  is  not  constant  along  the  ship's  length, 
only  the  wider  portions  of  the  ship  will  contribute  to  roll 
excitations.   It  is  easy  to. calculate  that  any  section  with 
less  than  60%  of  the  maximum  beam  will  have  a  negligible  con- 
tribution.  Therefore,  the  ship  can  be  considered  to  have  an 
effective  length  Le  which  should  be  used  in  the  formula 
L|cos(B)|  to  replace  the  waterline  length  L.   An  empirical 
formula,  which  is  a  conservative  estimate  based  upon  the  60% 
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max.  beam  criterion,  was  derived  graphically  using  a  Taylor 
Series  hullform  from  reference  9. 


1.0 

0.9 
0.8 
0.7 


L 


0.6 


L)  =  0.12  +  0.943  C 


P 


j i I I i ! I |_ 


0.6      0.7      0.8 


Figure  7.   Effective  Shiplength. 


The  region  of  validity  for  heading  angle  is  therefore  given 
by: 


Le|cos  (3)  |  <  A/2 


but  ...    X  =  2TTg/oo' 


Applying  this  standard  to  graphical  display  is  done  in 
Figure  8  below. 
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Figure  8.   Region  of  Validity  for  Heading  Angle. 

To  summarize,  it  is  obvious  that  the  range  of  parameter 
values  for  w  and  3  vary  significantly  with  ship  dimensions, 
ship  geometry,  and  sea  state.   It  is  important  in  each  case 
to  investigate  whether  the  extended  equations  are  still  valid 


7 .   Heading  Angle  Correction 

From  the  results  of  Section  4,  it  would  seem  desirable 
to  extend  the  range  of  heading  angles  further  by  some  approx- 
imate means.   An  easy  way  to  do  this  is  with  the  assumption 
that  the  static  righting  moment  is  uniformly  distributed 
along  the  effective  length  of  the  ship.   Then,  the  wave  exci- 
tation vanishes  anytime  that  Le | cos (3)  | /A  has  an  integer  value 
of  1  or  greater.   It  can  also  be  surmised  that  for  non- 
integer  values  of  Le | cos (8) | /A  greater  than  1,  the  part  of 
the  ship's  effective  length  between  successive  crests 
receives  a  net  zero  excitation  due  to  cancellation  effects. 
It  is  appropriate  to  include  this  effect  as  a  correction 
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factor  to  the  excitation  spectrum.   An  efficiency  parameter, 
e,  which  is  a  function  of  Le,  3,  and  X (w)  is  formed  by  put- 
ting the  above  arguments  into  mathematical  form. 


let:   X  =  Le|cos (3) |/X 


X  is  between 


0.0  0.5 

0.5  1.0 

1.0  1.5 

1.5  2.0 

2.0  2.5 

2.5  3.0 

3.0  3.5 

3.5  4.0 

4.0  4.5 


e[Le,  3,  X(w)] 


1.0 
2(1  -  X) 
2/3  (X-  1) 
2/3(2  -  X) 
2/5 (X  -  2) 
2/5(3  -  X ) 
2/7  (X  -  3) 
2/7  (4  -  X.) 


1.0 


Figure  9.   Heading  Angle  Correction 
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The  excitation  spectral  density  function  now  is 


.2  /r      2   .  2/n.   -2tbw2/g  „   ,  x 
\p   /6co  =  e   sin  (3)  e    u    3  Sf  f  (to) 


Since  an  actual  ship  will  probably  not  have  a  constant 
beam  along  its  effective  length,  K  will  not  be  uniformly 
distributed.   The  region  of  net  zero  excitation  will  not  be 
quite  equal  to  one  wavelength,  but  will  change  with  the 
relative  orientation  of  the  passing  wave.   Qualitatively, 
this  will  tend  to  introduce  higher  harmonics  into  the 
encounter  frequency  of  the  roll-inducing  wave  energy  than 
that  predicted  by  the  above  model.   However,  in  general  it  is 
safe  to  assume  that  the  cancellation  effect  does  occur,  even 
though  the  higher  harmonics  are  produced. 

Another  consideration  is  that  the  final  product  of 
interest,  roll  angle,  will  be  determined  from  integrating  the 
spectral  density  function,  and  taking  the  square  root  of  the 
result.   This  set  of  operations  makes  the  roll  angle  very 
insensitive  to  the  exact  form  of  the  heading  angle  correction 

Based  on  the  above  comments,  it  seemed  consistent  with  a 
conservative  approach  to  take  the  correction  factor  as  the 
line  connecting  the  peaks  of  the  function  shown  in  Figure  9. 

e  =  1.0,   X  <    0.5, 
e  =  1/2X,   otherwise. 
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8.   Spreading  Function 

A  spreading  function  was  used  to  convert  uni-directional 
long-crested  seas  (implied  by  S> J  into  a  more  realistic 
short  crested  sea  surface.   Following  reference  3,  the 
spreading  function  f  (y)  is  a  symetric  function  centered  on 
the  primary  direction  of  the  seas. 


f(y)  =  2/tt  cos  (y)  ,   -tt/2  <  y  <  tt/2 


S^2   f(y)S   (03)dy  =  1.0S   (co) 


The  effect  of  this  function  is  to  smooth  the  long 
crested  seas  prediction  for  roll  response  somewhat.   It  implies 
that  there  is  always  a  beam  component  of  wave  motion  influenc- 
ing the  ship  even  when  the  ship  is  headed  into  or  away  from  the 
predominant  direction  of  the  waves. 

9.   Sea  State  Data 

The  below  listed  values  were  taken  from  references  8,  11  and 
12,  and  are  useful  for  design  purposes.   The  peak  frequency 
of  the  surface  wave  elevation  spectrum  S>^  is  the  oo  given. 


p 

HV3 

Sea   State 

0.417 

30 

8 

0.455 

25 

7+ 

0.509 

20 

7 

0.589 

15 

6+ 

0.720 

10 

5 

1.02 

05 

3 

**** 

<05 

0-2 
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Probability  of  Occur- 
ence in  North  Atlantic 


0.01 
0.02 
0.03 
0.10 
0.35 
0.45 
0.04 
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IV.   Solution  for  the  Equations  of  Motion 

1.  Introduction 

Given  that  the  coefficients  are  all  known,  there  are  at 
least  three  practicable  methods  of  solution  available.   These 
are  analog,  numerical  integration,  and  complex-algebraic 
manipulation  of  the  linearized  equations.   Each  of  these 
methods  has  its  own  advantages  and  drawbacks  compared  to  the 
others.   The  latter  method  was  the  one  finally  adopted  because 
it  was  the  most  flexible  for  design  purposes  and  yields  the 
most  useful  information  for  effort  put  into  using  it.   How- 
ever, further  work  with  the  other  methods  may  be  warranted 
depending  upon  the  particular  anti-rolling  tank  problem  being 
investigated.   The  details  of  the  complex  algebraic  method 
are  given  in  this  chapter.   A  brief  overview  of  the  results 
of  the  investigations  into  the  other  two  methods  is  given  in 
Appendix  D. 

2.  Complex  Algebraic  Solution 

The  use  of  complex  algebra  along  with  the  spectral  den- 
sity technique  allows  for  a  closed  form  solution  of  the 
linearized  equations  of  motion.   The  basic  assumption  involved 
in  this  solution  bears  a  close  resemblance  to  what  may  be 
imagined  as  the  real  situation  of  the  response  of  a  ship  in 
waves.   It  should  be  recalled  that  gravity  waves  have  an 
advance  velocity  which  is  dependent  upon  frequency,  so  that 
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the  phasing  between  ocean  waves  is  essentially  random.   The 
waves  of  different  frequency  (or  wavelength)  travel  along  at 
different  velocities.   If  it  is  assumed  that  the  ocean  has  a 
constant  spectral  density  with  time,  then  the  ship's  response 
to  each  of  these  frequency  components  has  reached  a  steady 
state  value.   Random  phasing  of  these  components  accounts  for 
the  Rayleigh  type  distribution  of  peak  values  of  response, 
and  the  principal  of  linear  superposition  neatly  ties  the 
whole  system  of  equations  and  density  functions  together;  the 
ship's  response  is  the  sum  of  the  steady  state  responses  to 
each  of  the  frequency  components. 

A  mention  should  be  made  of  the  implications  of  the  term 
"linearity".   The  basic  mathematical  operation  for  each  com- 
ponent case  of  frequency  excitation  and  response  is  linear 
because  the  coefficients  have  been  determined.   However,  the 
expressions  for  the  so-called  linear  coefficients  which 
involve  w  ,  w  ,  X ,  phasing,  ship  dimensions  and  geometry, 
tank  saturation,  heading,  sea  state,  etc.  have  no  such 
restrictions.   They  are  required  only  to  be  analytical. 

The  standard  method  for  finding  the  steady  state  solu- 
tion to  a  set  of  ordinary  linear  differential  equations  is  to 
make  a  complex  transformation,  using  a  complex  algebraic 
operator  for  differentiation,  solve  the  complex  equations 
algebraically,  and  retrieve  the  information  of  response  mag- 
nitude and  phase  angle  from  the  complex  expression. 
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Let  ip  =  ^  e  ^  .   Then  the  steady  state  solution  must  be 


XZ       1-1-     C  A        A.  11JT  ljJT        * 

of   the    form:    $   =   <{>    e         ,    t\  =  T)   e    ^    ,    v* 


v*  e    ,  where 


if)    d>  -  n  ,  v    are  complex  constants.   The  operator 
Yo,      o   o   o         r  £ 

2      2 

D  =  d(  )/dx  may  be  replaced  by  D'=  iy,  D   =  -y  .   The  equa- 
tions of  motion  now  become  for  a  passive  tank: 


Cll(iy)  C12(iy)  C13(iy) 
C21(iy)  C22(iy)  C23(iy) 
C31(iy)       C32(iu)       C33(iy) 


o 


V 


o 


Lno 


1-X,  X    +  iX,  X 
by        bo 

X    +iX 
Y        o 


o 


c  (iy) 


v  *) 
o 


o 


1  -    X.  X      +    iX,X 
by  bo 

X      +    iX 
Y  o 


N(iy) 


determinant  [C(iy)] 


The  numerator,  N  (iy )  ,  is  found  by  replacing  the  first 

column  of  the  C  matrix  with  the  coefficient  matrix  for  \b    . 

ro 

Notice  that  the  C  matrix  terms  will  not  be  symmetric  because 
of  the  choice  of  v*  as  a  variable.   If  y*  were  chosen  instead 
of  v*,  the  symmetric  condition  for  an  energy  conservative 
system  would  have  been  met. 

Since  p  =  iy<|>  by  linear  theory,  the  roll  velocity 
response  may  be  obtained  by  y|~|  ,  and  likewise  the  roll 
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acceleration  is  y2!-^-!.   Tank  angle,  velocity  and  accelera- 
te 

tion  are  similarly  determined. 

The  final  step  in  obtaining  roll  response  is  to  multiply 
the  transfer  function  cj)2/^2  (or  §    2/ty2,    or  cf>  2/^2)  by  the 
excitation  spectral  density  function  i\)2/8ui   and  integrate  the 
result  over  the  entire  frequency  band.   The  spreading  func- 
tion f  (y)  may  be  inserted  in  this  step  depending  upon  whether 
long-crested  or  short-crested  sea  result  are  desired. 


/      CO 

<J)  =      /  /    ((f>2/ijj2 )  ^2/6w   doo  (long-crested    seas) 


V2  /— 

cf>  =      /      f  (y)   ///°°((})2/^2)^2/6lo   dwdy  (short- 

rms  -v/2  v    0  crested    seas) 


<J>1//3      =      2.0    <J) 

T  rms 
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V.   Anti-Rolling  Tank  Design 

1.   Introduction 

The  complex  algebraic  form  of  solution  of  the  non- 
dimensional  equations  of  motion  was  written  in  computer 
language  to  allow  the  numerical  work  to  be  done  by  a  digital 
computer.   Details  of  this  computer  program  are  given  in 

Appendix  C.   Essentially  there  are  two  separate  programs — 

1/3 
MAIN  1  which  maps  <f)  '       versus  u   and  K   for  one  case  of  head- 
ing angle  in  long-crested  seas  and  MAIN  2  which  calculates 
short-crested  sea  results  for  any  given  set  of  input  para- 
meters.  In  the  interest  of  numerical  efficiency  most  of  the 
iterative  steps  of  the  design  process  can  be  done  by  MAIN  1, 
looking  at  the  largest  component  of  roll  in  long-crested  seas 
which  makes  up  the  center  frequency  of  short-crested  seas 
results.   MAIN  2  can  be  used  later  in  the  final  stages  to 
give  short-crested  sea  roll  angle  predictions  based  upon  the 
results  of  using  MAIN  1. 

The  design  method  basically  entails  "bouncing"  back  and 
forth  between  MAIN  1  and  MAIN  2  looking  for  optimum  para- 
meters of  the  non-dimensional  equations.   The  steps  are  man- 
ually done,  meaning  that  each  step  represents  a  separately 
submitted  computer  run.   About  2  to  4  computer  runs  are 
required  for  a  final  design  for  each  different  anti-rolling 
tank. 
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The  passive  tank  design  method  in  section  2  suggests  how 
to  specify  various  parameters  of  the  equations  and  presents  a 
logical  sequence  to  arrive  at  an  approximate  tank  size. 
There  is  also  a  suggested  iterative  sequence  to  be  used  as  a 
design  aid.   One  of  the  fortunate  results  of  the  subsequent 
chapter  on  design  examples  is  the  collection  of  equations 
shown  on  Figure  12.   It  is  recommended  that  this  worksheet  be 
kept  as  a  tally  sheet  and  guide  when  performing  the  iterative 
steps. 

The  active  tank  design  method  given  in  section  3  is 
mostly  a  discussion  of  possible  changes  to  the  passive  tank 
design  which  are  compatible  with  the  computer  program.   Mak- 
ing these  changes  will  result  in  a  prediction  of  the  perfor- 
mance of  a  feasible  activated  anti-rolling  tank. 

2.   Passive  Tanks 

The  response  of  the  ship-passive  tank  system  is  modeled 
by  the  non-dimensional  equations  developed  previously  and 

setting  E  (t)  =0. 

3   o 

Tank  location  may  be  specified  or  a  small  number  of  dif- 
ferent locations  may  be  investigated.  The  reason  for  limit- 
ing tank  location  to  the  absolute  minimum  is  entirely  practi- 
cal. The  vertical  space  of  the  ship  is  divided  by  structure 
into  a  number  of  decks  and  for  access  and  strength  considera- 
tions the  cross  duct  of  the  U  tube  tank  should  run  across  (or 
under)  a  deck.   Also  there  may  only  be  a  single  or  a  few 
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regions  within  the  ship  large  enough  to  accept  the  tanks.   As 
a  consequence  of  the  non-dimensional  equations,  there  is  a 
curious  coupling  between  tank  vertical  placement,  size,  and 
the  coefficients  of  the  equations.   Since  the  ship's  vertical 
center  of  gravity  is  shifted  with  any  vertical  shift  of  the 
tank,  the  ship's  GM  is  reduced  whenever  the  tank  is  moved 
higher  in  the  ship.   Since  K   is  proportional  to  GM,  and  oj 
is  proportional  the  square  root  of  K  ,  and  because  all  of  the 
non-dimensional  coefficients  depend  on  K,  or  w  ,  each  verti- 
cal  tank  location  must  be  treated  as  a  separate  problem.   The 
design  method  which  follows  is  for  each  separate  case  of  ver- 
tical tank  location. 

Once  tank  location  is  specified,  the  parameters  which 

are  still  free  to  be  chosen  are:  A^ ,  y, ,  y  ,_ ,  K^ .   We  know 

t    t    st    t 

that  K^  has  a  minimum  value  of  about  -0.05,  and  that  y,  will 
t  t 

have  a  value  close  to  1.0  because  the  vibration  absorber 
effect  of  the  tank  is  maximized  by  tuning  the  tank  to 
y.  =  w./w   =  1.0.   These  four  dimensionless  parameters  are 
dependent  upon  the  real  ship  and  tank  parameters  as  follows: 


*4.  =  X.  [K,  (GM)  ,pA  R2,m] 
t     t   cf>        o   ' 


yt  =  ytta)s(V  '  R'S'  (A0/Al'Lo/R)1 


yst  =  yst[ws(K  } 'R'S"(zt/R*Lo/R)] 

2 

Kt  =  K  [pA  R  ,  geometry  of  tank,  A  /A,] 
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If  the  ship  changes  its  loading  condition  during 
expected  operations  such  that  its  displacement  is  appreciably 
changed,  then  all  four  parameters  have  a  A  dependency,  also. 

Another  important  parameter  for  design  is  ty   max.   Tank 
capacity  ty   max  is  defined  as  the  maximum  heel  angle  that  the 
tank  can  impart  to  the  ship  by  maximum  movement  of  the  tank 
fluid  off  to  one  side.   This  also  represents  the  maximum 
waveslope  angle  which  the  tank  is  capable  of  stabilizing.   To 
form  an  opinion  of  a  required  ip  max,  use  must  be  made  of 
empirical  formulas.   Hagen  mentions  the  use  of: 


\b  =    0  .  36/Log.,  .  (A)  radians 

rmax        '       ^10 


A  formula  suggested  by  center  of  buoyancy  correction  to  the 
spectral  density  function  which  takes  the  equivalent  of  a  6' 
surface  waveslope  at  the  center  of  buoyancy  for  the  ship's 
natural  resonant  frequency  is: 


*    =  0.10  e"(tbWs2/g)  radians 
max 


The  above  formulas  give  an  idea  of  the  approximate  capacity 
which  a  tank  must  have  in  order  to  be  effective,  and  serve  as 
a  feasibility  check  before  proceeding  further  with  a  particu- 
lar design. 

2 

The  parameter  X      =  (-2pgA  R  /-mgGM)  has  certain  char- 
acteristics which  allow  it  to  be  specified  by  the  designer. 
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2 

Rearranging  the  formula  to  read  A   =  [2A  R  /(displaced  volume 

of  ship)]/(GM),  notice  that  the  numerator  of  this  ratio  is 
the  standard  free  surface  correction  for  a  virtual  rise  in 
the  ship's  center  of  gravity  due  to  the  fluid  being  unconfined 
in  the  n.  direction.   Therefore,  A,  is  the  percentage  loss  in 
static  stability  of  the  ship  due  to  the  tank  fluid  free  sur- 
face.  There  will  be  a  maximum  allowable  A   which  ensures 
safe  operation  of  the  ship  under  all  loading  conditions, 
including  damaged  conditions.   Since  a  larger  A   is  associ- 
ated with  more  tank  stabilizing  effectiveness  (tj>    -A^L  /R)  , 

3  rmax   to    ' 

the  choice  of  A   is  fixed  at  the  maximum  value  which  allows 
safe  operation  of  the  ship. 

An  important  conclusion  to  reach  from  the  above  observa- 
tions is  that  the  ship  hull  form  should  be  designed  with  the 
notion  that  an  anti-rolling  tank  is  to  be  installed.   This 
means  that  space,  weight,  and  static  stability  will  be  such 
that  the  ship  will  be  comfortable  with  a  tank  of  the  capacity 
indicated  by  the  empirical  formulas.   The  alternative  choice 
of  trying  to  fit  an  anti-rolling  tank  into  an  existing  ship 
design  which  may  not  have  all  the  necessary  reserve  space, 
weight,  or  static  stability  most  likely  will  result  in  an 
ineffective  stabilizing  system. 

Now  that  A,  is  known,  and  \b  has  been  determined 

t  max 

approximately,  an  assumption  for  a  y   value  is  all  that  is 

required  to  calculate  tank  dimensions  (A  ,  R,  L  ,  z^)  the 

o      o    t 

tank  parameter  y   (L  /R,z  /R),  and  approximate  tank  weight 
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[2  pA  R(L0/r  +  A./A  )].   A  first  guess  at  y   is  made  by  using 
Chadwick & Klotter ' s  equation,  based  upon  the  linear  equations, 
neglecting  sway,  from  reference  4. 


w  =  u  -  Atd/yEt2  -  i/Mst4>r1/2 


It  is  now  possible  to  calculate  all  the  tank  dimensions 
and  all  the  coefficients  of  the  equations  of  motion  with  suf- 
ficient accuracy  to  consider  them  constant,  except  for  duct 

area  An  and  tank  parameters  y,  and  K , . 
1  t       t 

The  remainder  of  the  design  problem  may  be  stated  as 
follows:  Find  the  best  combination  of  y   and  K  which  mini- 
mizes the  rolling  motion  predicted  by  the  non-dimensional 
equations  of  motion.   This  two  parameter  search  problem  may 
be  solved  by  noticing  that  y   must  be  about  1.0,  and  that  two 
digit  accuracy  is  the  best  to  be  expected  from  the  equations 
of  motion.   Therefore,  instead  of  a  sophisticated  search 
technique  it  should  suffice  to  map  roll  angle  versus  y   and 
K, .   This  may  be  further  simplified  by  sub-optimizing  each 

y   value  for  K  .   This  approach  is  justified  by  gentle  cur- 
l/3 
vature  of  the  function  cj>    (y.,K  ),  in  the  region  of  interest, 

as  found  in  all  sample  problems  tried.   The  remainder  of  this 

section  will  investigate  numerical  methods  to  deal  with  this 

search  problem. 

The  basic  strategy  for  solving  this  problem  is  to  first 

fix  upon  some  criterion  of  roll  performance.   Based  upon  the  sea 
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state  data  already  given,  the  author  decided  to  minimize  roll 
angle  at  the  worst  heading  at  operational  speed  in  sea  state 
6  (15  feet  significant  wave  height)  was  a  reasonable  criter- 
ion.  Other  criteria  could  include  roll  velocity  or  roll 

acceleration  for  other  environmental  conditions. 

1/3 
In  order  to  compute  <f>     for  given  conditions,  a  long 

1/3 

sequence  of  operations  is  accomplished--^     is  twice  the 

2   2         2 
square  root  of  the  integral  of  <J)  /y      times  m  /Sco.   This 

sequence  makes  vivid  the  reason  why  the  complex  algebraic 

method  of  finding  cj>/ip  is  superior  to  numerical  integration. 

The  latter  would  require  another  integration  of  many  steps,  a 

very  expensive  approach  since  the  results  must  in  turn  be 

iterated  many  times  for  a  good  design. 

The  above  discussion  assumed  that  long  crested  seas  were 

involved.   For  short  crested  seas  another  integration,  this 

time  over  every  possible  heading  angle,  must  be  performed  for 

the  spreading  function.   For  this  reason,  long  crested  seas 

1/3 
are  assumed  for  mapping  (J)    [y  ,  K  (best)  ]  .   Short  crested 

seas  are  only  looked  at  later  as  an  output. 

1/3 
In  mapping  cj)    [u  ,  K  (best) ]  ,  it  was  found  that  the 

function  displays  very  gradual  curvature  as  was  to  be 

expected  for  an  integral  function.   Therefore,  only  a  few 

points  are  investigated  during  each  computer  run,  because 

interpolation  is  predictable.   A  second  order  Legendre 

polynomial  curve  is  used  in  an  interpolation  scheme  to  find 

a  local  minimum  or  maximum.   Manual  inspection  of  the  computer 
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results  will  show  whether  further  mapping  is  necessary. 
Warning:  the  interpolation  function  likes  to  extrapolate  to 
find  minimums  or  maximums  equally  well,  but  such  results  are 

invalid  by  the  statement  of  the  problem. 

1/3 
Once  <J)    [u.  ,  K  (best)  ]  has  been  mapped  and  values  of 

ll  and  K^  for  a  minimum  roll  angle  have  been  decided,  short- 
t      t 

crested  sea  results  may  be  obtained.   The  problem  then  shifts 
to  chosing  which  of  the  alternative  tank  locations  is  the 
best.   Economic  as  well  as  performance  factors  will  enter 
into  the  decision  of  this  question.   For  example,  a  tank 
which  is  in  a  cargo  hold,  presumably  valuable  real  estate, 
would  probably  not  be  chosen  over  one  which  was  in  a  less 
costly  part  of  the  ship,  such  as  atop  an  outside  deck,  were 
the  two  tanks  to  perform  nearly  as  well. 

A  detailed  listing  of  the  author's  steps  in  passive  tank 
design  are  reported  as  an  aid  in  improving  upon  them.   A 
later  section  will  include  a  completely  worked  example  of  the 
method. 


I .   Problem  setup 

a)  Assume  a  range  of  a.  small  number  of  vertical  tank 
locations . 

b)  Calculate  ty  requirement. 

max    ^ 

c)  Specify  a  X^ 

t 

d)  Assume  an  initial  u,  value  from  Chadwick  and  Klotter's 
equation. 
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e)  Determine  Lo,  A  ,  R,  z  ,  and  approximate  A,. 

f )  Design  the  tank  except  for  A, . 

g)  Calculate  the  hydrodynamic  coefficients  and  other 
computer  inputs. 

II .  First  Iteration 

a)  Run  MAIN  2  program;  use  the  long  crested  wave  results 
to  estimate  the  worst  speed  and  heading  angle  com- 
bination.  Use  y   assumed,  K   =  -0.8 

1/3 

b)  Map  (J)    [y  ,  K  (best)  ]  using  MAIN  1.   Use  y   assumed, 

K   initial  =  -0.4,  Ay   =  0.5,  AK   =  -1.0. 

1/3 

c)  Interpolate  to  find  cf>     minimum.   Check  computer 

results  to  ensure  that  a  true  minimum  occurred. 

III.  Second  Iteration 

a)  Run  MAIN  2  again  to  see  if  the  worst  speed  and  head- 
ing angle  combination  has  changed.   Also,  roll  per- 
formance may  have  degenerated  at  other  headings  and 
speeds  significantly. 

b)  Continue  iterations  if  change  in  speed  or  heading 
angle  warrant  it. 

IV.  Alternative  Tank  Choice 

a)  Based  upon  the  operator's  or  owner's  criteria,  select 
the  best  tank  location. 
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V.   Verification 

As  a  final  step,  short  crested  sea  results,  MAIN  2, 
should  be  run  for  the  entire  expected  range  of  speed  and 
sea  conditions.   This  is  to  make  sure  that  the  solution 
didn't  fix  too  shortsightedly  upon  some  assumed  worst 
condition  and  lessen  it  without  overlooking  its  negative 
effects  on  some  peculiar  or  unforseen  combination. 
Another  verification  which  should  be  performed  during 
this  and  all  other  steps  as  well  is  to  check  the  range  of 
parameters,  the  range  of  the  solution  variables  to  ensure 
that  their  coefficients  remain  linear  (or  else  adjust  the 
coefficients  to  their  equivalent  linear  values)  and  that 
the  \i      and  K   values  being  predicted  are  realizeable. 

3.   Activitated  Anti-Rolling  Tanks 

An  activated  anti-rolling  tank  is  one  which  has  an  input 
power  source  in  addition  to  ship  induced  excitations.   A 
pressure  source  is  used  as  the  type  of  power  input  in  this 
discussion,  as  opposed  to  a  flow  source.   It  is  desired  to 
regulate  the  pressure  source  in  such  a  way  as  to  reduce  the 
rolling  motion  of  the  ship. 

The  choice  of  a  pressure  source  instead  of  a  flow  source 
was  not  entirely  arbitrary.   A  flow  source  would  have  to  be 
some  sort  of  a  positive  displacement  pump  capable  of  handling 
large  volume  rates  of  fluid  transfer,  approximately  500  tons 
per  minute  as  an  order  of  magnitude.   The  required  pump  size 
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for  this  volumetric  capacity  would  make  it  economically 
unfeasible.   Another  requirement  for  the  power  source  that  is 
imposed  by  economics  and  amply  illustrated  in  other  papers, 
e.g.  reference  18,  is  that  it  operate  at  a  nearly  steady 
power  level.   The  peak  power  output  for  a  pressure  source,  as 
found  by  standard  control  theory  methods  using  linear  feed- 
back, is  too  high  to  be  practical.   There  is  too  much  wasted 
power  capacity  which  is  only  used  at  peak  power  demand.   How- 
ever, the  pressure  pump  must  be  built  large  enough  to  be  able 
to  supply  the  predicted  peak  demand  if  the  linear  equations 
are  to  be  valid.   Once  further  economic  factor  to  be  con- 
sidered is  the  competition,  the  alternative  methods  of  roll 
damping.   With  regards  to  these,  the  gyroscopically  controlled 
anti-rolling  fin  installation  is  the  most  serious  candidate. 
In  practice,  a  standard,  well  designed  anti-rolling  fin 
installation  has  been  found  to  consume  about  one  or  two  per- 
cent of  the  installed  hoursepower  of  a  ship,  including  pro- 
pulsion, due  to  the  added  effects  of  the  fin  area  increasing 
the  ship's  total  hull  resistance  to  passage  through  the  water 
and  the  power  needed  by  the  fin  positioning  motors.   Thus, 
economics  and  practical  engineering  considerations  have  dic- 
tated a  rather  narrow  range  of  possible  control  effectors  for 
anti-rolling  tanks — a  pressure  type  source  of  power  of  nearly 
constant  output  which  has  a  power  rating  of  about  one  percent 
of  the  ship's  installed  power. 
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A  further  consideration  of  the  power  source  is  reversi- 
bility.  A  level  power  output  does  not  imply  that  the  direc- 
tion of  pressure  or  flow  will  be  constant.   The  pressure 
supply  will  obviously  have  to  reverse  itself  as  the  ship  rolls 
back  and  forth.   This  aspect  of  the  pressure  source,  quanti- 
fied by  a  characteristic  time  lag  to  change  direction,  will 
have  to  be  reasonably  estimated  and  included  in  the  analysis, 
and  must  also  be  reasonable  in  the  sense  that  a  smaller  time 
lag  probably  represents  a  more  expensive  piece  of  equipment. 

It  is  helpful  to  first  look  at  the  orders  of  magnitude 
of  the  parameters  needed  to  describe  a  control  effector. 
Consider  a  ship  of  20,000  tons  displacement  with  a  speed  of 
20  knots.   Assuming  an  overall  lift-to-drag  ratio  of  150,  the 
installed  horsepower  would  be  about  18,400  Hp.   This  would 
indicate  from  the  preceeding  arguments  that  an  economically 
feasible  anti-rolling  tank  can  be  no  greater  than  184  Hp. 
Compare  this  to  the  energy  flow  of  the  ship  in  its  roll 
coordinate,  assuming  a  GM  of  6 ' ,  a  5°  roll,  and  a  roll  period 
of  12  seconds. 

20,000  tons  (2240  lb/ton) (6 ') (sin  5°) (1/2) • 
(1/12  sec)  (Hp-sec/550  ft-lb)   =   1775  Hp 

This  simple  estimate  is  enough  to  show  that  we  cannot 
expect  the  control  effector  to  be  able  to  counter  the  effects 
of  wave  excitations.   However,  it  may  be  used  to  counter  the 
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resonance  phenomenon  associated  with  rolling  frequencies  near 
the  natural  resonance  peaks  of  the  ship-tank  system. 

Since  the  passive  anti-rolling  tank  is  capable  of  devel- 
oping nearly  as  much  energy  flow  as  the  ship,  until  the 
saturation  angle  of  the  tank  is  reached,  it  seems  justifiable 
to  further  simplify  the  control  effector  equation  to 


T   p(t)  +  p(t)  =  ±P 

p  r      r        c 


P   denotes  the  steady  pump  pressure  developed  across  the  tank 
duct,  and  T   is  the  reversing  lag  time  for  the  pump.   The 
equation  P   =  f  {$  ,  cf>  ,  n ,  r\)    implies  now  that  the  sign  for  P   is 
chosen  according  to  some  criterion  contained  in  the  state 
variables  (\>,<\>fr\,    and  r\.         Consider  the  same  ship  with  an 
anti-rolling  tank  of  200  tons  weight  of  fluid.   The  average 
fluid  volume  flow  would  be: 


0.707(200  tons)  (35  ft3/ton  H20)  (1/12  sec)  =  142  ft3/sec 


For  an  average  184  Hp,  the  mean  pressure  is  184  Hp  (550  ft-lb/ 
Hp-sec)  (1/142  ft3/sec)  =  245  lb/ft2  =  1.7  psi.    This  pressure 
changes  the  fluid  level  between  the  tanks  by  3.8  ft.   Assum- 
ing a  50  percent  pump  efficiency,  the  pressure  differential 
would  be  0.85  psi  with  a  fluid  level  difference  of  1.9  ft. 
This  corresponds  to  a  tank  angle  of  about  two  degrees.   The 
static  ship  roll  angle  for  still  water  and  constant  pump 
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pressure  would  be  A   times  two  degrees,  amounting  to  only  a 
fraction  of  a  degree.   Thus  we  can  now  state  that  the  addi- 
tion of  such  a  control  effector  would  have  a  qualitatively 
negligible  effect  on  the  equations  of  motion,  except  near  the 
resonant  frequencies.   This  indicates  that  some  form  of 
active  damping  is  the  goal  of  the  control  effector  switching 
criterion . 

The  problem  of  adding  a  control  effector  now  becomes  one 
of  how  to  add  the  active  damping  to  the  linearized  equations 
of  motion.   The  tank  equation  is  the  only  one  of  the  set 
which  contains  p(t)  for  an  active  tank.   Substituting  P   = 
g  (J)  +  g  <j)  +  g  n  +  g.r)  into  the  tank  equation,  neglecting  for 

the  moment  T   (consider  T   small) ,  and  combining  terms  in  the 
P  P 

linearized,  non-dimensional  equations  of  motion,  the  struc- 
ture of  the  active  tank  equations  is  revealed. 


2     2    Wsg2           gl          D2 
t     Kst     2pgR  2pgR   r    r-j  v* 

yyt 

+  A,  [D  /\i.       +  ( —  -  ~ — )D  +  (1  -  ~ =r)  ]  n   =   0. 

t      t      u     2pgR'  2pgR 


We  would  expect  that  the  effects  of  g..  and  g^  to  be 
relatively  small  even  if  we  took  them  at  their  maximum  pos- 
sible values  (  to  develop  184  Hp)  because  of  the  small  static 
values  of  tank  and  roll  angle  the  pump  is  able  to  develop. 
Hence,  making  g,  =  g3  =  0  is  approximately  correct  for  any 
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economical  tank.   The  choice  of  g„  and  g.  to  use  as  a  switch- 
ing criterion  is  obvious.   The  tank  already  has  linear  damp- 
ing which  may  be  varied  as  required  by  adding  valves  or 
baffles.   Negative  tank  damping  'is  undesirable  because  it 
would  increase  roll  response  at  resonance  frequencies.   The 
computer  program  results  from  MAIN  1  confirms  this  hypothesis 
Adding  a  damping  term  to  the  ship  roll-tank  coupling  terms, 
however,  would  have  a  desirable  effect  at  any  frequency.   An 
additional  tank  moment  caused  by  the  pressure  differential 
would  be  induced  opposing  the  roll  velocity  of  the  ship. 
Thus,  the  controller  effects  may  be  included  in  the  linear- 
ized equations  of  motion  for  the  type  of  controller  described 
by  substituting  this  reduced  tank  equation  into  the  non- 
dimensional  equations  of  motion. 

2     2   Wsg2  XtD 

yyt 

D2     KtD 
+  AM— 2  +  -^-   +  l]n   =   0 
t  Pt     Ut 

The  value  of  g2  is  determined  by  power  considerations. 
The  energy  dissipation,  the  product  of  assumed  pump  pressure 
and  tank  volume  flow,  can  be  calculated  and  compared  to  the 
desired  pump  power  in  an  iterative  fashion  until  the  two 
balance.   The  effect  of  pump  lag  time  may  be  included  by 
deducting  the  fraction  of  the  period  of  roll  that  was  spent 
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in  changing  direction.   This  factor  is  about  2tt/oo  -  T   for 

T   less  than  2tt/w  and  zero  for  T   greater  than  2tt/w. 
P  P 

The  above  active  tank  equations  are  only  useful  for  pre- 
dicting ship  motions  in  a  seaway  which  overpowers  the  passive 
aspects  of  the  anti-rolling  tank.   Fortunately,  this  is  the 
case  when  an  active  tank  would  most  likely  be  used.   The 
additional  damping  term  contained  in  g~  and  the  need  to  cal- 
culate mean  pump  power  are  small  modifications  to  the  com- 
puter program  for  the  equations  of  motion.   Otherwise,  the 
same  design  method  as  that  for  the  case  of  passive  anti- 
rolling  may  be  used. 

The  type  of  hardware  which  would  give  the  desired  pump 
characteristics  is  a  waterjet  type  pump  situated  in  the  tank 
cross  duct.   The  jet  velocity  would  be  much  higher  than  the 
mean  fluid  velocity,  giving  rise  to  a  steady  pressure  inde- 
pendent of  average  tank  fluid  velocity.   The  presure  rise  is 
caused  by  the  cross  duct  acting  as  an  inefficient  diffuser 
and  the  pressure  reverse  may  be  accomplished  by  switching  the 
waterjet  outlet  from  a  nozzle  in  one  direction  of  flow  to  an 
opposite  nozzle.   The  actual  pump  motor  would  then  not  have 
to  reverse  at  all,  and  the  lag  time  would  be  associated  with 
the  time  required  to  switch  nozzles. 
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VI.   Design  Examples 

In  order  to  validate  the  computer  program,  verify  the  design 
steps  postulated,  and  investigate  the  possible  range  of  the  solu- 
tions, it  was  decided  to  design  a  series  of  antirolling  tanks 
for  a  sample  family  of  ships.   The  ships  were  chosen  to  be  able 
to  illustrate  as  many  aspects  and  implications  of  the  equations 
of  motion  as  possible  within  time  and  effort  constraints. 

Three  ships  were  finally  settled  upon  to  show  basically  the 
effects  of  ship  size,  degree  of  sway  coupling,  and  tank  size. 
They  were  a  3600-ton  small  naval  escort  type,  a  14,500-ton  medium 
sized  naval  vessel  which  resembles  a  small  aircraft  carrier 
or  amphibious  landing  ship,  and  a  20,400-ton  commercial  stores 
carrier.   The  three  ships  were  assumed  to  be  existing  designs 
to  which  an  antirolling  passive  tank  was  to  be  added.   The  two 
naval  vessels  were  possessed  of  a  large  degree  of  sway  coupling 
introduced  by  the  parameter  Ab  =  BG/B  because  they  were  ships 
with  high  centers  of  gravity  relative  to  their  other  dimensions. 
The  commerical  stores  carrier  had  only  about  half  as  much  sway 
coupling,  being  more  typical  of  a  cargo  vessel  which  does  not 
possess  the  high  topside  payload  weights  of  a  naval  vessel. 

1.   Ship  Size 

The  principal  dimensions  of  the  three  vessels  were  as  follows: 

Ship  No.  1    Ship  No.  2    Ship  No.  3 
A  3600         14,500        20,400 

L  408  585  500 
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— G 


Flight  Deck 


Upper  Gallery  Deck 


LOY/er   Gfll.lery   Deck 


Hangar  Deck 


1 


'igure  10B  Midship  Section  Ship  No.  2,  With  Tank 
Locations 
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Ship  No.  1     Ship  No.  2  Ship  No.  3 
B                 45           76.67        71.4 
T               14.5           21.33        28.6 
Cb           0.4733         0.5304        0.70 


2.   Tank  Size 

Based  upon  the  empirica] formula  of  Section  V.2,  an  approxi- 
mate required  value  for  tank  capacity,  ^max.,  was  determined. 
The  vertical  layout  of  the  ships,  by  decks,  is  shown  in  Figure  10 
Tank  dimensions  and  parameters  were  calculated  using  the  collec- 
tion of  equations  shown  on  figure  12  by  an  iterative  process. 
The  A.  value  for  all  tanks  was  maximized  in  the  following  way: 
An  assumed  A,  was  used  for  the  low  tank  and  applied  to  the  value 
of  metacentric  height,  GM,  for  the  low  tank.   The  residual  stabi- 
lity of  the  reduced  metacentric  height  given  by  GM  (1  -  A  ) 
was  kept  as  a  required  value  to  be  achieved  by  the  other  tank 
locations.   Since  the  higher  tank  locations  resulted  in  a  higher 
center  of  gravity  for  the  ship  and  therefore  lower  GM,  their 
A,  values  had  to  be  necessarily  lower  in  order  to  achieve  the 
required  residual  stability  GM  (1  -  A  ) . 

Ship  No.l    Ship  No.  2   Ship  No.  3 


¥  max.    =    .36/logio    A 

5.8° 

5.0° 

4.8° 

.„                       ,    -tbu)2/q 
V  max.    =    .le           s'  y 

5.4° 

5.4° 

5.2° 

GM    (1    -    A,., 

2.69 

5.17 

5.25 

Of  the  three  deck  locations  possible  for  tanks  on  Ship  No.  1, 
it  was  found  that  the  tank  placed  on  the  main  deck  had  too  small 
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a  capacity  ¥  max.  and  ship-tank  coupling  A   to  be  practical. 

Only  the  lower  tank  and  the  center  tank  were  further  investigated, 

On  Ship  No.  2  there  were  only  two  tank  locations  possible, 
assuming  that  keeping  the  center  hangar  deck  area  clear  was 
imperative.   The  higher  tank  was  slightly  cramped  in  order  to 
fit  it  between  decks,  while  the  lower  tank  should  probably  be 
expanded  to  meet  the  hangar  deck  bottom  to  eliminate  the  small 
wasted  space  between  the  tank  top  and  the  hangar  deck. 

On  Ship  No.  3,  the  presence  of  vertical  cargo  access  hatches 
and  cargo  storage  limited  again  the  tank  location  to  two  possi- 
bilities:  a  low  tank  with  the  cross  duct  underneath  the  cargo 
floor,  or  to  a  second  deck  high  tank  of  small  dimensions  to 
accommodate  the  cross  duct  between  cargo  hatches. 

For  the  naval  vessels,  Ships  No.  1  and  2,  the  tank  fluid 

was  assumed  to  be  fuel  oil  such  as  diesel  oil  or  JP-5  (  P  =  1.69 

3 
slugs/ft.  )  because  such  vessels  carry  a  generous  reserve  fuel 

supply  for  extended  range  purposes.   In  effect,  this  made  the 
passive  antirolling  tanks  a  reserve  fuel  oil  tank  which  is  pos- 
sible to  ballast  with  seawater  in  an  emergency.   For  the  com- 
mercial vessel,  Ship  No.  3,  the  tank  weight  was  too  large  to 

consider  it  to  be  carrying  expensive  fuel  oil,  so  seawater  (  P  = 

3 
1.99  slugs/ft.  )  was  assumed  to  be  the  tank  fluid.   This  re- 
duced the  tank  volume  somewhat,  compared  with  a  fuel  filled 
tank,  but  the  tank  weights  of  between  1-3  percent  of  the 
ship's  total  displacement  was  a  factor  against  their  being  a 
practical  commercial  installation. 
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3.  Operational  Criteria 

The  choice  of  speed  and  sea  state  was  also  made  to  illus- 
trate the  range  of  solutions  of  the  equations  of  motion.   The 

1/3 

smallest  ship  was  examined  in  a  sea  state  5  (H    =  10  ft.) 

at  its  operational  speed  of  25  knots,  a  typical  requirement 

for  an  escort  ship  of  its  class.    Ship  No.  2  was  looked  at 

1/3 
in  sea  state  6  (H  '   =15  ft.)  at  a  lower  speed  of  20  knots 

because  it  can  be  expected  to  behave  better  than  its  lighter 

counterpart,  and  has  a  lower  designed  speed.   Ship  No.  3  was 

examined  at  its  economical  cruising  speed  of  16  knots  in  a 

1/3 

storm  of  sea  state  7  (20  ft.  H  '     ) .   This  is  because  the  tank 

should  be  able  to  mitigate  motion  effects  on  the  cargo  in  ex- 
treme conditions. 

4 .  Initial  Program  Values 

An  initial  estimate  of  the  worst  heading  angle,  the  point 
for  optimizing  roll  angle  versus  y,  and  K  ,  was  made  by  taking 
the  point  where  the  peak  of  the  excitation  spectral  frequency 
distribution  corresponded  to  the  higher  frequency  peak  of  the 
roll  response,  assumed  about  130  percent  of  go  .   Later  work 
showed  that  this  was  not  consistently  a  very  good  first  guess, 
and  that  a  better  first  estimate  could  be  gotten  by  running 
Main  2  program  with  Kt  =  -0.8  and  taking  the  worst  heading  angle 
directly  from  the  results. 

The  initial  values  for  heading  angle  and  other  parameters 
for  the  first  iteration  were  as  follows  for  each  of  the  ship- 
tank  combinations. 


71 


Ship 

No.  1 

Ship 

No.  2 

Shi 

p  No.  3 

low 

center 

low 

high 

low 

high 

BG 

9.75 

9.93 

16.5 

17.3 

8.37 

8.69 

GM 

3.36 

3.18 

7.39 

6.67 

7.00 

6.68 

fcb 

5.67 

5.67 

8.24 

8.24 

13.43 

13.43 

At 

.20 

.155 

.30 

.22 

.25 

0.075 

yt 

1.01 

1.01 

1.02 

.97 

1.01 

1.00 

2 

yst 

4.11 

4.10 

4.72 

-24.61 

4.25 

-88.3 

Kp 

-.221 

-.221 

-2.92 

-2.92 

-3.17 

-3.17  (X  107) 

|Yym| 

.67 

.67 

.53 

.53 

.88 

.88 

|K^/Ix| 

.26 

.26 

.41 

.41 

.15 

.15 

Dl 

-8250 

-8250 

-30200 

-30200 

0 

0 

D2 

-21450  • 

-21450 

-78500 

-78500 

-96790 

-96790 

D3 

-10300  • 

-10300 

-37700 

-37700 

-26800 

-26800 

Fl 

0 

0 

0 

0 

0 

0 

F2 

.684 

.684 

.684 

.684 

.735 

.735 

F3 

1.05 

1.05 

1.05 

1.05 

1.13 

1.13 

3 

90° 

9  0° 

118° 

118° 

50° 

50° 

u 

43 

43 

34 

34 

27 

27 

Hl/3 

10 

10 

15 

15 

20 

20 

KG 

18.58 

18.76 

30.4 

31.1 

23.54 

23.86 

K<J> 

-2.71   ■ 

-2.57 

-24.0 

-21.66 

-31.99 

-30.52  (X  107) 

ws 

.5413 

.5271 

.4456 

.4233 

.505 

.493 

Ao 

124.5 

91.3 

607 

357 

611 

214 

Ao/Aj^ 

5.0 

5.26 

4.61 

4.85 

3.4 

3.83 

R 

20 

20 

33 

35 

32 

32 

Lo/R 

.498 

.538 

.303 

.286 

.349 

.313 
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SHIP:                                                                              TANK: 

\b          -  At   Lo/R 
max 

iteration 

iteration 

iteration 

iteratior 

K.    =  -mg   x   GM 
9 

-2pgAoR2 

u  - 

Lo  =    (R/At)    i|i 

max 

A  K, 

A               t    * 
o   "'  -2pgR2 

/                    1 

^st2          yst4 

/     -\ 

"S        /m(.38B2(l+|^|) 

lx. 

0)      =    u    to 
t          t  s 

1 

A      /  A                                                         t        / 

AQ/A     -          2              ~^/R 
tot   R/g 

03           -      /£                      X 

st        /      R  2Lo/R  +  Z    /R 

y       =  a)     /w 
St           st      s 

_   2pgA  R 

t          99/n      (L   /R  +  A./A   ) 
224U           o               1      o 

KG1    -    (AKG  +   6W5z)/A 

GM'    =   KG  +  GM  -   KG' 

Figure  12.   Tank  Sizing  Equations  -  Worksheet 


Ship 

No.  1 

low 

center 

Weight 

84 

64.6 

zt/R 

.343 

-.339 

wst 

1. 

.097 

1.067 

76 

Ship  No.  2       Ship  No.  3 

low     high  low     high 
506       299     611       163  (tons) 
.436      -.78     .23      -.67 

.968     i2.1  1.041    i4.75 


5.   Summary  of  Results 

1/3 
The  results  of  calculations  for  tj>  '   as  a  function  of  y.  at 

K.  (best)  are  shown  on  figure  13  for  the  first  three  iterative 
steps  of  looking  at  an  assumed  point  in  the  sea   state  -  heading 
angle  space.   A  discussion  of  results  for  each  of  the  three 
ships  follows  this.   It  should  be  repeated  that  the  intent  was 
not  to  design  a  particular  antirolling  tank  or  to  apply  any 
arbitrary  criteria,  but  instead  to  get  as  broad  a  look  as  pos- 
sible at  the  range  and  implications  of  the  equations  of  motion, 
and  the  solution  method. 

Because  the  center  tank  of  Ship  No.  1  did  not  seem  to  offer 
any  better  performance  than  the  low  tank,  it  was  not  carried 
beyond  the  first  iteration.   This  low  tank  in  Ship  No.  1  became 

the  most  difficult  of  the  tanks  to  design  because  the  plot  of 

1/3 

cf>  '  *  versus  y   at  K.  (best)  exhibits  more  curvature  and  other 

surprises  than  the  other  tanks.   At  each  iteration,  the  sub- 
optimized  parameters  which  worked  well  at  that  particular  head- 
ing angle  were  poor  at  different  headings.   Finally,  y.  was 
chosen  as  1.01  by  inspection  of  all  the  data,  and  K.  was  deter- 
mined by  running  Main  2  program  for  short  crested  seas  through 
Kfc  =  -0.4,  -0.8,  and  -1.0.   The  combination  y   =  1.01,  K   =  -0.8, 
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which  was  the  first  iteration  guess,  was  just  about  optimum. 
Overall  performance  of  the  antirolling  tank  compared  to  the 
ship  without  a  tank  was  not  spectacular,  but  because  the  tank 
uses  reserve  fuel  oil  as  a  working  fluid  and  is  partially  sit- 
uated in  normal  tankage  space,  it  was  a  reasonably  successful 
design. 

On  Ship  No.  2  the  low  and  high  tanks  were  both  carried 
through  to  a  final  iterative  step.   The  initial  guess  of  118 

as  a  worst  heading  angle  delayed  the  solution  somewhat  because 

1/3 
at  that  heading  the  characteristic  of  $   '       did  not  flatten  out 

as  expected  and  also  118   was  nowhere  close  to  the  real  worst 
heading.   Both  high  and  low  tanks  performed  well  and  it  is 
interesting  to  note  that  in  spite  of  the  much  smaller  size, 
weight,  and  A   of  the  high  tank  it  performed  as  well  as  the  low 
tank.   This  was  probably  because  the  parameter  u  . 2  for  the 
high  tank  had  a  negative  sign,  indicating  that  roll  acceleration 
pushed  the  tank  fluid  in  the  direction  to  give  a  positive  tank 
righting  moment.   For  y   2  with  a  positive  sign,  the  case  for 
any  tank  below  the  ship  center  of  gravity,  the  opposite  reac- 
tion of  tank  fluid  to  roll  acceleration  occurs. 

On  Ship  No.  3  only  the  low  tank  was  investigated  fully, 
but  it  can  be  noticed  at  the  50   heading  angle  that  the  high 
tank's  performance  is  nearly  as  good  as  that  of  the  larger  low 
tank.   Again,  this  was  due  to  the  negative  sign  on  jj   2. 

A  relative  comparison  of  the  predicted  solutions  of  the  most 
important  tank  design  parameters  is  in  order. 


83 

Ship  No.  1    Ship  No.  2    Ship  No.  2    Ship  No.  3 
low  tank     low  tank     high  tank     low  tank 


*P  max 

5.7° 

5.2° 

3.2° 

5.0° 

Xt 

0.20 

0.30 

0.19 

0.25 

*t 

1.01 

1.02 

1.20 

1.01 

yst 

2.03 

2.72 

i5.30 

2.06 

Weight/1 

00A 

2.33 

3.49 

1.57 

3.00 

Kt 

-0.8 

-1.0 

-1.0 

-.4 

£ 

0.73 

0.26 

0.25 

0.26 

1/3  1/3 

A  performance  parameter  e  (e  =(j)  '      max  (tank) /(J)  '   max  (no  tank)) 

was  added  at  the  bottom  of  the  list,  but  it  is  difficult  to  attach 
much  relative  meaning  to  the  parameter  because  all  of  the  ships 
had  different  operational  criteria.   Even  if  the  criteria  were 
the  same,  the  size  differential  in  the  ship's  displacements  would 
have  made  comparisons  meaningless  except  in  the  case  of  the  tanks 
on  the  same  ship.   The  tank  capacity,  fmax,  did  not  enter  into 
the  optimization  because  all  of  the  significant  tank  angles  cal- 
culated were  less  than  their  saturation  angles,  the  arctan  (L  /R) . 
Of  the  parameters  given  above,  the  first  four,  ^max,  y,  ,  A.  ,  and 
u  .  are  the  most  fundamentally  important  because  they  determine 
the  antirolling  tank  size,  dimensions,  weight,  placement  in  the 
ship  vertically,  and  they  cannot  be  changed  once  the  installation 
is  completed.   K,  however  may  be  readily  adjusted  to  suit  the 
particular  ship  and  the  other  two  parameters  relate  to  performance 
or  economic  matters. 
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6.   Analysis  of  Results 

a.   Tank  Capacity.   The  question  of  how  well  the  tank  capa- 
city was  utilized  may  be  answered  by  inspection  of  the  results 

of  significant  tank  angle  at  the  worst  heading  in  short-crested 

1/3 
seas.   The  predicted  tank  angle,  n  '  ,  may  be  compared  to  the 

tank  saturation  angle,  arctan  (Lo/R) ,    to  give  a  measure  of  how 

much  capacity  was  used  at  the  design  conditions.   Two  tank  angle 

1/3 

statistics  were  looked  at.   n    is  a  good  engineering  estimate 

of  the  point  at  which  saturation  effects  must  be  taken  into 
account  in  calculations.   n  '     was  arbitrarily  taken  as  an 
estimate  of  when  tank  saturation  effects  might  just  start  to  be 

noticed  by  the  crew  of  the  vessel  because  of  the  slapping  of  the 

1/3 
tank  fluid  against  the  top  of  the  tank.   n    is  an  output  of 

the  MAIN  2  program  and  n      is  found  by  the  statistics  of  the 

1/3 
random  process  to  be  3.34/2.0  0  times  n    . 


n1/3 
nl/100 


arctan  (Lo/R) 

,1/3/aI 
1/100 


n^/arc tan  (Lo/R) 


Ship  No.  1 

Ship  No.  2 

Ship  No.  2 

Ship  No.  3 

low  tank 

low  tank 

high  tank 

low  tank 

9.81° 

7.07° 

7.24° 

13.48° 

16.38° 

11.81° 

12.09° 

22.51° 

26.47° 

16.86° 

15.96° 

19.24° 

0.47 

0.53 

0.58 

0.89 

R)   0.62 

0.70 

0.76 

1.17 

An  immediate  observation  is  that  the  low  tank  on  Ship  No.  3 
is  the  only  one  which  is  making  good  use  of  its  tank  capacity  at 
design  conditions.   The  other  designs  are  all  wasting  space  by 
building  too  large  a  tank  for  the  requirements  stated. 
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b.  A  .   The  ship-tank  coupling  parameter  did  not  by  itself 
seem  to  influence  results  much  except  as  it  contributed  to  the 
tank  capacity.  Ship  size  and  tank  placement  had  a  more  definite 
effect. 

c.  y, .   The  tank  tuning  parameter  had  a  secondary  effect 
on  the  results  because  a  higher  u   corresponding  to  a  higher  co 
is  realized  by  decreasing  the  area  of  the  cross  duct,  Al .   The 
lighter  cross  duct  weight  allowed  for  more  of  the  fluid  to  be 
in  the  wing  tanks.   Consequently,  a  change  in  tank  tuning  to  a 
higher  tank  frequency  made  for  a  tank  of  slightly  greater  capacity, 

d.  u   .   The  ship-tank  decoupling  frequency  parameter  was 
important  for  the  reasons  previously  discussed  concerning  it. 

A  more  negative  value  of  1/y  .  2,  which  is  accomplished  by  moving 
the  tank  higher  vertically,  is  always  associated  with  better  tank 
performance  for  a  tank  of  a  given  A   and  capacity  fmax. 

e.  Weight  Fraction.   This  parameter  can  be  directly  related 
to  the  cost  of  the  tank.   A  complementary  parameter,  the  space 
fraction,  was  not  given  because  the  data  was  unavailable  for 
useful  ship  volume. 

f.  K.  .   The  tank  damping  parameter  falls  within  the  narrow 
range  of  about  -0.4  to  -1.4  for  all  successful  designs.   This 
statement  is  based  on  analysis  of  much  more  data  than  is  presented 
in  the  results  section.   A  convincing  reason  for  the  truth  of  this 
statement  can  be  found  by  considering  the  critical  damping  for 

a  single  degree  of  freedom  system.   For  a  single  degree  of  freedom 
system  critical  damping  implies  that  there  is  no  overshoot  and 
no  oscillatory  response  to  a  step  input.   Such  critical  damping 
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has  the  K   value  of  -2.0.   Therefore  the  ratio  of  tank  damping 
to  critical  damping  falling  in  the  range  of  0.2  to  0.7  is  rea- 
sonable if  the  tank  is  to  act  as  a  dynamic  vibration  absorber.   A 
K,  prediction  of  greater  than  -2.0  tells  us  that  the  given  tank 
parameters  are  not  consistent  with  a  stabilizing  tank.   An 
example  of  K.  greater  than  -2.0  is  shown  on  figure  15  for  the 
low  tank  on  Ship  no.  1. 

g.   Tank  Placement.   There  is  another  aspect  of  vertical 
tank  placement  in  addition  to  comments  already  given  about  how 
it  affects  X.  ,  y  .  and  GM.   A  special  case  of  tank  operation  is 
whenever  the  duct  is  closed  for  some  reason  while  the  tank  is 
full,  for  example  for  emergency  maintenance.   On  high  tanks 
especially  there  is  a  large  GM  loss  due  to  the  KG  rise  associated 
with  the  large  tank  fluid  weight  high  in  the  ship.   Figure  16A 
gives  for  an  example  the  high  tank  on  Ship  No.  2  to  show  the 
undesirable  effects  of  this  operational  condition. 

h.   Ship  Size.   It  seems  obvious  from  the  results  to  state 
that  a  lighter  displacement  ship  such  as  Ship  No.  1  apparently 
is  not  as  well  suited  a  candidate  for  an  antirolling  tank,  instal- 
lation as  a  larger  ship.   Whether  this  is  due  to  the  center  of 
buoyancy  being  less  submerged  or  to  the  inertia  of  the  ship 
allowing  for  a  higher  natural  frequency  remains  to  be  answered. 

7 .   General  Conclusions 

Although  none  of  the  design  examples  is  complete  in  itself  - 
primarily  because  a  range  of  speeds  and  sea  states  was  not  looked 
at  for  each  ship  -  it  was  possible  to  draw  a  number  of  general 
conclusions  from  them. 


87 

First  of  all,  the  equations  of  motion  seem  to  be  valid. 
Uncoupling  sway  by  setting  Xb  =  0  and  1/u  , 2  =  0  gives  the  clas- 
sical  fourth  order  differential  equation  for  a  coupled  ship-tank 
system.   The  assumed  excitations  are  the  primary  departure  from 
the  latter,  along  with  the  sway  coordinate,  and  hence  introduce 
the  greatest  possible  source  of  error  because  the  reactive  force 
and  moment  terms  were  also  taken  from  classical  theory.   A  lot 
of  numerical  accuracy  hinges  upon  the  two  critical  assumptions 
of  uncoupling  of  roll  and  sway  and  application  of  the  excita- 
tions at  the  center  of  (static)  buoyancy.   Model  testing  of  some 
sort  is  needed  to  prove  or  disprove  these  conjectures. 

The  general  design  method  seems  to  be  flexible  and  effi- 
cient within  the  range  of  interest.   This  was  due  to  the  veri- 

1/3 
fication  that  <$>  versus  y   and  K ,  ,  a  'drawable'  three  dimen- 

sional function,  is  a  relatively  flat  function.   This  was  due 
to  its  being  an  integral  function. 

The  solution  of  the  fifth  order  differential  equations  of 
motion  by  the  design  method  predicts  tank  parameter  values  which 
are  quite  similar  and  in  some  cases  identical  to  the  classical 
fourth  order  approach.   The  roll  responses  of  the  two  approaches 
is  not  significantly  different.   However,  a  more  refined  solu- 
tion is  possible  by  the  methods  of  this  paper.   For  example,  the 
high  tank  of  Ship  No.  2  with  u.  =  1.2  had  a  lower  weight  than 
that  for  the  classically  designed  tank  of  y,  =  0.97  due  to  the 
smaller  cross  duct  associated  with  a  higher  tank  frequency. 
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VII.   Conclusions  and  Recommendations 

1.   Conclusions 

The  design  method  proposed  in  this  paper  seems  to  offer  an 
attractive,  efficient  and  flexible  method  for  the  design  of  pas- 
sive antirolling  tanks.   The  analysis  for  active  antirolling  tanks 
led  instead  to  a  method  for  a  feasibility  study  for  the  use  of 
active  tanks.   The  active  tank  design  method  predicts  perform- 
ance of  an  economical  active  tank  and  will  select  optimum  tank 
parameters  for  it.   However,  the  design  of  the  pump  and  the  pump 
controller  was  left  as  a  separate  step.   The  active  tank  analysis 
also  demonstrated  why  active  tanks  are  such  unpopular  creatures. 
The  small  difference  in  tank  angle  which  the  pump  is  able  to 
achieve  and  the  small  amount  of  active  damping  which  it  adds  to 
the  equations  of  motion  are  not  nearly  enough  to  overcome  the 
chief  cause  of  tank  efficiency  suffering  in  high  sea  states  — 
the  saturation  of  the  tank  which  places  an  upper  limit  on  the 
tank  righting  moment. 

The  design  method  for  passive  antirolling  tanks  is  attrac- 
tive because  it  uses  a  good  representation  of  the  real  ocean 
wave  process  in  predicting  roll  angle.   Since  there  is  much 
available  information  on  ocean  conditions  for  specific  ship- 
ping routes  or  operational  areas,  this  information  may  be  sub- 
stituted for  the  assumed  Pierson-Moskowitz  function.   The  design 
method  is  efficient  because  the  equations  of  motion  lend  them- 
selves to  a  complex  algebraic  form  of  solution,  thereby  resulting 
in  a  high  numerical  efficiency  for   the  computer  program.   Because 
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of  this  a  large  range  of  solutions  may  be  investigated  at  little 
cost.   Finally,  the  design  method  is  flexible  in  that  tank  de- 
sign may  be  approached  from  a  wide  variety  of  given  requirements. 
For  example,  meeting  some  criteria  for  operational  use,  meeting 
space  or  weight  constraints,  or  any  combination  of  these  may  be 
imposed  as  design  requirements.   Also,  the  equations  for  the 
hydrodynamic  and  tank  coefficients  may  be  changed  as  necessary 
in  the  main  computer  program  without  altering  the  subroutine 
which  calculates  roll  response. 

The  prediction  of  roll  angles  seems  reasonably  well  in  agree- 
ment with  all  of  the  references,  both  in  magnitude  and  in  dis- 
tribution over  heading  angles.   The  short-crested  sea  predictions 
should  be  more  accurate  than  those  for  long-crested  or  regular 
waves  because  linearity  is  better  preserved.   Only  small  compo- 
nents of  roll  for  each  heading  angle  are  added  together  by  super- 
position, including  the  worst  heading  at  resonance  peaks. 

A  general  conclusion  concerning  design  parameters  is  that 
a  reasonable  passive  tank  design  for  preliminary  design  purposes 
may  be  estimated  for  a  given  X,   ,  or  A   =  0.20,  u  ,  from  Appendix  B, 
u.  from  the  Chadwick-Klotter  equation  on  page  53,  K.  =  -0.80, 
and  the  tank  dimensions  from  Appendix  B. 

2.   Recommendations 

It  was  decided  to  first  investigate  and  validate  the  equa- 
tions of  motion  by  taking  the  simplest  case  of  a  passive  tank 
and  linearized  coefficients.   As  it  became  apparent  that  the  re- 
sults of  this  approach  were  promising,  a  design  method  for  passive 
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tanks  was  devised.   A  thorough  checking  of  the  results  by  model 
or  full  scale  testing  seems  to  be  in  order  as  the  next  step. 

The  numerical  integration  technique  seems  to  offer  the  best 
hope  for  modelling  a  first  order,  or  higher  order  control  effector 
in  conjunction  with  the  fifth  order  ship  and  tank  equations. 
However,  as  previously  mentioned  this  technique  is  inefficient 
in  that  there  are  so  many  integration  steps  to  be  performed  in 
order  to  do  an  optimization,  and  several  orders  of  magnitude 
more  integration  steps  to  be  able  to  look  at  short  crested  seas 
over  a  range  of  sea  states  and  ship  speeds.   In  short,  the 
numerical  analysis  technique  has  a  certain  merit  for  analysis 
but  is  a  slumsy  synthesis  tool.   Now  that  the  basic  equations 
for  the  reactive  forces  and  moments  have  been  validated,  there 
may  be  a  way  to  overcome  this  difficulty.   This  is  to  combine 
the  spectral  density  function  and  the  spreading  function  before- 
hand to  make  again  an  irregular,  short-crested  sea  by  taking 
all  the  harmonic  frequencies  and  superposing  them.   A  five  or 
ten  minute  real  time  segment  of  excitations  could  be  made  and 
stored  as  library  data.   This  data  could  be  manipulated  for 
speed,  heading  angle  and  sea  state  to  form  a  simulation  of  the 
excitations.   This  could  be  used  to  drive  the  equations  of  motion 
with  the  effector  included.   Data  for  events  such  as  significant 
roll  angle  could  be  obtained  from  the  statistics  of  the  results. 
It  would  be  interesting  to  see  how  well  the  results  of  such  a 
simulation  compare  with  the  complex  algebraic  solution.   There 
is  room  for  question  as  to  whether  the  simulation  would  be  as 
accurate  as  the  complex  algebraic  solution  because  of  the  reduced 
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form  some  of  the  linear  coefficients  would  assume.   For  example, 
the  sway  damping  term  must  lose  its  frequency  dependence  in 
this  type  of  approach  because  the  excitation  would  be  a  series 
of  absolute  quantities  derived  from  a  frequency  spectrum  and 
would  lack  the  frequency  information  necessary  for  calculation 
purposes. 

It  was  noticed  when  comparing  the  results  of  a  passive  tank 
design  with  the  design  solution  contained  in  reference  13  that 
there  was  very  close  agreement  in  overall  results  in  short- 
crested  seas  except  for  the  case  of  the  region  of  heading  angles 
near  0   and  near  180  ;  that  is,  head  seas  and  stern  seas.   The 
author  believes  that  the  discrepancy  most  likely  could  be  due 
to  the  simplistic  approach  taken  towards  obtaining  a  heading 
angle  correction  in  Section  III . 7.   This  author  assumed  a  constant 
distribution  of  wave  excitations  over  an  effective  shiplength. 
It  is  recommended  that  other  distribution  functions  more  closely 
resembling  that  which  is  actually  found  aboard  the  ship  be  tested 
to  see  whether  a  closer  agreement  could  be  realized.   The  impli- 
cations of  using  a  second  or  third  order  parabolic  distribution 
would  be  worth  examining.   Before  doing  this,  however,  it  would 
be  wise  to  carry  through  on  the  first  recommendation  given  in 
this  section  and  check  on  which  of  the  predictions  of  roll  angle 
are  closer  to  model  or  full  scale  experimental  tests,  those  of 
this  author  or  those  of  reference  13. 

A  final  recommendation  is  that  some  of  the  steps  for  passive 
tank  design  could  be  better  automated.   For  example,  the  inte- 
grals S1  and  S"  could  be  accurately  handled  for  any  situation 
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by  the  proper  numerical  program.   Also,  the  burden  of  solving 
the  tank  sizing  equations  on  Figure  14  which  need  an  iterative 
convergent  solution  should  be  removed  from  the  tank  designer. 
This  would  be  no  simple  task  of  programming,  however,  because 
of  the  great  number  of  combinations  of  constraints  such  as  co 
or  a)    being  limited  by  tank  dimensional  constraints. 


93 


References 


1.  Abkowitz,  M.  A.,  Stability  and  Motion  Control  of  Ocean 

Vehicles .   MIT  Department  of  Ocean  Engineering, 
Cambridge,  Mass.,  1969 

2.  Bell,  J.  and  W.  P.  Walker,  "Activated  and  Passive  Controlled 

Fluid  Tank  System  for  Ship  Stabilization".   The 
Society  of  Naval  Architects  and  Marine  Engineers  Paper 
No.  4,  1966 

3.  Comstock,  J.  P.,  ed. ,  Principles  of  Naval  Architecture, 

Chapter  IX,  The  Society  of  Naval  Architects  and 
Marine  Engineers,  New  York,  1967 

4.  Chadwick,  J.  H.  and  K.  Klotter,  "On  the  Dynamics  of  Anti- 

Rolling  Tanks",  Supplement  to  Division  of  Engineering 
Mechanics,  Stanford  University,  Technical  Report  No. 
15,  N6-ONR-251-2,  February  16,  1953. 

5.  Chadwick,  J.  H. ,  "On  the  Stabilization  of  Roll",  Journal 

of  the  Society  of  Naval  Architects  and  Marine  Engi- 
neers, vol.  63,  1955 

6.  Chu,  Dalzell,  and  Modisette,  "Theoretical  and  Experimen- 

tal Study  of  Ship-Roll  Stabilization  Tanks",  Journal 
of  Ship  Research,  September,  1968 

7.  Dalzell,  J.  F. ,  "A  Note  on  the  Distribution  of  Maxima  of 

Ship  Rolling",  Journal  of  Ship  Reserach,  December, 
1973 

8.  Gale,  P.,  "Naval  Surface  Ship  Design;  Selection  of  Stabili- 

zation System  Type  and  Capacity" ,  Summer  Lecture 
Series,  MIT  Department  of  Ocean  Engineering,  August, 
1974 

9.  Gertler,  M. ,  A  Reanalysis  of  The  Original  Test  Data  for 

the  Taylor  Standard  Series,  David  Taylor  Moden  Basin 
Report  No.  806,  Navy  Department,  March  1954. 

10.  Hagen,  G.  R. ,  Feasibility  Studies  of  the  Roll  Stabilization 

of  the  USS  BOSTON  (CAG-1) ,  David  Taylor  Model  Basin 
Report  No.  950,  Navy  Department,  September,  1955. 

11.  Kehoe,  J.  W. ,  "Destroyer  Seakeeping,  U.S.  and  U.S.S.R.", 

Naval  Engineer's  Journal,  December,  1973 


94 


12.  Loukakis,  T.  A.,  Computer  Aided  Prediction  of  Seakeeping 

Performance  in  Ship  Design,  MIT  Department  of  Ocean 
Engineering,  Report  70-3,  August,  1970 

13.  Miller,  E.  R.  and  Lain,  H.  Comparative  Performance  of 

Roll  Stabilization  Systems  for  the  PF,  Fly  Dronautics, 
Inc.  Technical  Report  No.  7103-3,  March,  1972 

14.  Neuman,  J.  N.  ,  Marine  Hydrohynamics ,  Second  Revised  Edi- 

tion, MIT  Department  of  Ocean  Engineering,  1972 

15.  Reuter,  W. ,  "Naval  Surface  Ship  Design;  Lecture  47-48", 

Summer  Lecture  Series,  MIT  Department  of  Ocean 
Engineering,  August,  1974 

16.  Stigter,  C. ,  "The  Performance  of  U-Tanks  as  a  Passive 

Anti-Rolling  Device",  Nederlands  Scheepsstudiecen- 
trum  TNO,  Report  81S,  Delft,  October,  1964 

17.  Van  Den  Bunt,  J.  D.,  "The  Design  of  U-Tanks  for  Roll 

Damping  of  Ships",  Nederlands  Scheepsstudiecentrum 
TNO  Report  124S,  Delft,  October,  1969 

18.  Webster,  W.  C. ,  "Analysis  of  the  Control  of  Activated 

Antiroll  Tanks",  The  Society  of  Naval  Architects 
and  Marine  Engineers  Paper  No.  9,  1967 


95 


Appendix  A  Anti-Rolling  Tank  Geometry  Factors 


It  was  decided  to  construct  a  parametric  family  of  Anti- 
Rolling  tanks  using  the  simple  tank  geometry  of  Figure  A-l.   In 
practice,  a  much  more  complicated  geometry  could  conceiveably 
exist  due  to  the  need  to  use  the  ship's  bulkheads,  sides,  and 
decks  as  part  of  the  tank  structure.   Also,  depending  upon  the 
tank  damping  required,  it  may  be  necessary  to  fair  in  sharp 
corners  or  other  flow  obstructions.   It  is  expected,  however, 
that  the  assumed  geometry  will  be  close  enough  to  any  required 
geometry  in  any  specific  installation  so  that  equivalent  dimen- 
sions and  equivalent  integrals  (S1  and  S")  can  be  quickly  cal- 
culated and  overall  dimensions  easily  estimated. 

The  basic  geometry  is  shown  below. 

I 
G 


2Lo+%  h 


Figure  A-l.   Anti-Rolling  Tank 
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The  tank  form  factors  are  S1  and  S" . 
S'  =    j(Vo/a)dl 

S"  =    J(d/R)dl 

Performing  the  required  integrations  by  parts  over  each 

of  the  regions  of  constant  sectional  area,  A  ,  A, ,  and  A  ,  and 
_>  '   o    1       o 

dropping  the  nearly  negligible  higher  order  terms  gives  the 
desired  results. 

S'  =  2(LQ  +  (R  -  bt)AQ/A1)  =  2R(LQ/R  +  Aq/A;l) 

S"  =  2R(2L  /R  +  z./R) 
o      t 

The  two  tank  natural  frequencies  are  to   and  to   . 


wt  =V2g/S' 


u)st  =V2g/S' 


0.) 


t  "V^"  •  -\//(Lo/R  +  Ao/Al) 


wst  v^yv^Lo^,  2t/R) 

The  above  steps  have  reduced  the  problem  of  selecting  a 
tank  size  to  selecting  a  mix  of  the  tank  parameters  R,  (A  /A, ) , 
(L  /R) ,  and  (z ,/R) .   It  may  be  seen  that  the  other  tank  dimen- 
sions are  known  once  these  are  established.   The  range  of 
parameters  for  which  this  analysis  is  valid  is  approximately: 

0.2  <  LQ/R  <  1.0 

-  1  <  zt/R  <  +  2 

R  <  B/2 

Notice  that  for  z.  <-2L   the  value  of  to  .  will  be  imaginary 

to  st  3     J 

Later  work  with  the  equations  show  that  this  indicates  a  90 
phase  shift  of  the  motion  induced  in  the  tank  by  the  ship  roll 
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motion,  but  does  not  imply  any  dramatic  change  in  the  equations 
of  motion. 

Figure  A-2  summarizes  the  relationship  between  the  tank 
frequencies  and  the  tank  parameters.   For  most  combinations  of 
these  parameters,  and  in  the  range  of  validity  previously  men- 
tioned, Figure  A-2  is  sufficiently  accurate  for  preliminary 
tank  design  purposes. 

A  convenient  way  to  determine  tank  size  is  by  fluid  weight. 

Tank  weight  (fluid)  =  2pA  R(L  /R  +  A,/A  ) 
^  o    o       1   o 

An  effectiveness  paraneterof  tank  stabilizing  potential  is 
the  ratio  of  tank  righting  moment  to  wave  upsetting  moment  per 
unit  of  tank  or  waveslope  angle.   This  ratio,  A. ,  is  also  equal 
to  the  percentage  reduction  in  static  GM  that  the  tank  induces 
on  the  ship. 

At  =  2pgAoR2/K<t>  =   1 2A   R   /35AGM   (A  in  long  tons,  *=  p/pH  0) 
Another  measure  of  tank  stabilizing  potential  is  the  maxi- 
mum waveslope  that  the  tank  could  statically  stabilize.   Since 
the  maximum  tank  angle  possible  is  about  L  /R  (radians) ,  the 


static  tank  capacity  i|>  max  is  given  by 


\\)  =  AtL  /R 

rmax      o' 


\p  =  2KA  L  R/35AGM 

max      o  o 
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Appendix  B.   Coefficient  Estimates 

1.   Hydro-dynamic  Coefficients 

Figures  B-l  and  B-2  display  the  hydrodynamic  coefficients 
for  ellipsoidal  bodies  of  revolution,  as  theoretically  pre- 
dicted by  potential  flow  theory.   The  strategy  in  using  these 
figures  is  to  assume  a  body  with  the  same  dimensions  (L,B,T,) 
as  the  ship.   In  cases  where  the  mass  of  the  body  of  revolution 
is  quite  a  bit  different  than  that  of  the  actual  ship,  correction 
can  be  made  by  taking  the  average  value  of  the  coefficient  for 
a  body  of  the  same  dimensions  (L,B,T)  and  for  a  body  of  the 
same  mass  as  the  ship  with  similar  geometry  (L/B,  B/T,  L/T) . 

A.  Kp  and  Yv 

These  are  both  constants.   For  order  of  magnitude 

estimates,  Yv  =  -0.9m,  and  K.  =  -0.25  I  . 

P  x 

B.  K   and  Y 
_P v 

K  ,  given  as  B44  on  the  figure,  is  assumed  to  be  a 
constant  mean  value  for  reasons  previously  discussed.   B~~  on 
the  figure  corresponds  to  Y    (w,U)  by  the  notation  of  this 
paper.   It  is  assumed  to  be  well  enough  represented  as  a  se- 
quence of  linear  segments.   Figure  B-3  graphically  depicts  this 
linear  segmented  approximation.   The  coordinates  (F-,  D.)  are 
used  as  computer  program  inputs.   The  function  is  assumed 
constant  after  point  3;  although  this  is  certainly  grossly 
inaccurate,  the  roll  response  operator  for  the  normal  range  of 
ships  is  so  small  in  this  region  that  any  assumption  about  the 
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behavior  of  the  actual  function  will  have  essentially  no  effect 
on  the  roll  angle  solution. 

The  sway  damping  component  Y  ~ (U)  is  estimated  by  the 

methods  of  reference  (3)  to  be: 

2  2 

-TTpT  U  _<  Yv2(U)  £  -(Trp/2)  T  U 

The  lesser  absolute  value  is  associated  with  a  bare  hull 
only,  and  the  factor  of  two  is  added  for  a  ship  with  a  large 
amount  of  rudder  surface  area  and  other  stabilizing  underwater 
sail  areas  such  as  skegsf  keels,  etc.   A  particular  ship  may 
be  somewhere  between.   For  a  normal  amount  of  rudder  area  it 
will  probably  be  close  to  the  larger  absolute  value. 

C.   K£ 

K<|>   may  be    found   directly    from   the   formula   K<J>   =   -mg   GM. 

2.   Ship  Parameters 

A.  t,  ,  the  center  of  buoyancy 

The  center  of  buoyancy  is  accurately  estimated  by 

using  Morrish's  Formula  for  KB. 

t,  =  T  -  KB 
b 

KB  =  1/3(5/2  T  -  35A/Awp) 
Awp  =  Cwp  •  L  •  B 
Cwp  =  0.237  +  0.833  Cp   reference  14. 

B.  GM  and  G. 

The  vertical  center  of  gravity  and  GM,  if  not  known, 
may  be  estimated  together  by  the  relations: 
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KG  =  KB  +  BM  -  GM 

BM  =  Cit  L(B)3/35A 

Cit  =  -0.013  +  0.10  Cp   reference  14 

GM  =  0.08  B  (see  reference  3  for  a  list  of  repre- 
sentative GMs) 

Notice  that  if  one  or  the  other  of  G  or  GM  are  known,  the 

unknown  quantity  is  very  accurately  determined.   Otherwise, 

both  quantities  suffer  from  being  subject  to  'guesstimates'. 

C.   Ix 

The  transverse  moment  of  inertia  of  the  mass  of  the 

2 

ship  may  be  closely  estimated  by:  k   =  0.38  B;  Ix  =  m  (k  ) 
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Appendix  C.   Computer  Program 

1.   General  Information 

A  program  which  provides  a  numerical  solution  to  the 
equations  of  motion  was  written  in  standard  Fortran  language. 
It  is  compatible  with  any  Fortran  capable  system  as  long  as 
the  input  device  number  KI  and  the  output  device  number  KO  are 
correctly  specified  in  the  program.   A  plotting  subroutine 
PICTR  was  used  in  the  program;  if  this  library  subroutine  is 
unavailable,  replace  the  CALL  PICTR  cards  with  dummy  CONTINUE 
cards. 

There  are  two  basic  programs,  Main  1  and  Main  2, 
which  are  separable  decks  and  need  not  be  run  concurrently. 
A  single  subroutine  RESP   calculates  the  roll  angle  per  unit 
waveslope,  the  tank  angle  per  unit  waveslope,  and  the  phase 
angle  between  the  tank  and  roll.   A  library  subroutine,  PICTR, 
is  called  to  give  a  pictorial  display,  but  because  the  results 
are  printed  out  anyway,  it  is  unessential  to  the  program. 

Main  1  section  of  the  program  computes  optimum  tank 
damping  and  tuning  for  a  given  set  of  conditions. 

Main  2  section  of  the  program  computes  roll  response 
versus  heading  angle  in  short  crested  seas  for  a  given  ship, 
tank,  and  sea  state. 

Because  there  are  many  untried  variations  to  the 
program  it  is  advised  that  anyone  using  it  should  proceed  with 
care  and  check  results  at  each  step  of  the  computations. 
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2.   Computer  Program  Symbols 

The  following  symbol  list  cross  references  computer 
variables  with  the  symbols  used  elsewhere  in  this  text,  where 
indicated. 

MAIN  PROGRAM 


IX 

Ix 

LB 

Ab 

KFI 

K(J> 

KPD 

Kp 

KY 

K 

y 

KS 

K 
s 

LY 

Ay 

KP 

K 
P 

H3 

Hl/3 

LE 

Le 

C 

Spreading  Func 

Coefficients 

XLAB 

Axis  Labels 

TITLE 

Title  Tables 

XMT 

lit 

XXMT 

Dummy  ut 

F13 

*1/3 

RAO 

Response  Ampl. 

Operator 

S  or  SK  \\)    2/6  u> 

AB        Dummy  holder 


O 


U) 


E 

Tank  Angle  n 

TANG 

n1/3 

KO 

Output  device  number 

KI 

Input  device  number 

BG 

BG 

GM 

GM 

CB 

Cb 

XL 

L 

B 

B 

T 

T 

TB 

fcb 

XYVD 

| Yv/m | 

XLT 

At 

YV2 

Y 
V 

U 

U 

XKT1 

K.  initial  value 

DKTl 

K.  increment 

XMST2 

ust 

MAIN  PROGRAM  (continued) 
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DMT  ut  increment 

Wl  co  initial  value 

DW  co  increment,  Aco 

XKPD  |Kp/Ix| 

BETA  3 

Dl  -  D3   Sway  Damping  Y  , 
values 

Fl  -  F3   co  values  correspon- 
ding to  Di 


YV 
YVD 

Gl 


Yv 


Y. 
v 


U 


G2 

/2 

WS 

ws 

WE 

we 

W 

0)0 

XI 

X 

EP 

Heading  angle  cor- 

rection, 1/2  X 

CTB 

-2tbco2//g 
e       ^ 

ETA  or  E 

n 

YV1 

Y 
V 

N 

Number  of  increments 

of  Aco 

LT 

At 

L0 

Ao 

LY 

Ay 

LB 

Ab 

M 

y 

M2 

2 

y 

MT2 

yt 

MST2 

Vis 

KY 

Ky 

KS 

K 
s 

KT 

K^ 

SUBROUTINE  RESP 

W 

coo 

WS 

cos 

Gl 

fl 

G2 

a 

TOCBL 

Phase  angle 

Aij 

Matrix  elements , 

real 

Bij 

Matrix  elements, 

imaginary 

Ci 

Intermediate  matrix 

steps 

RX 

♦/* 

RT 

n/* 
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3.   Computer  Program  Input 

The  following  input  data  are  required  by  the  program. 

Ship  Dimensions 

L  Ship  Length 

B  Ship  Beam 

T  Ship  Draft 

Cb         Block  Coefficient  (Immersed  volume/LBT) 

BG         Distance  from  center  of  gravity  to  center  of  buoy- 
ancy, +  in  z  direction 

GM         Metacentric  Height 

t.  Draft  to  the  center  of  buoyancy 

Le         Effective  shiplength  (between  60%  beam  sections) 

Tank  Parameters 
ut         Tank  tuning  (cot/ws) 

Apt        Tank  tuning  increment  for  search  routine 
Kt  Tank  damping 

Ar.         Tank  damping  increment  for  search  routine 

2 

ust        Ship-tank  decoupling  frequency;  primarily  a  function 

of  vertical  placement 

At         Ship-tank  coupling  parameter  --  %  loss  in  static 
GM  of  ship  due  to  tank  free  surface  effect 

Equation  Coefficients 
Roll  damping 

Added  mass  coefficient  in  sway. 
Added  inertia  coefficient  in  roll. 
Sway  damping  ordinates  for  Y  ,  (u),U) 
Sway  damping  corresponding  to  D. 


K 
P 

lVml 

|Kj/Ix| 

Dl    -    D3 

Fl    -    F3 
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col 

Aw 


H 

3 

u 


1/3 


Excitation  Parameters 

Initial  value  (lower  cut-off  frequency,  about  0.2) 

Increment  of  u>  (notice,  30  increments  are  used  un- 
less N  is  changed  in  the  program;  the  30  increments 
plus  a),  will  determine  the  upper  cut-off  freq.) 

Significant  waveheight  of  seas. 

Heading  angle  of  ship  with  respect  to  the  seas 
(180   for  head  seas) 

Ship  speed  in  the  x-direction 


Program  Control  Variables 

There  are  three  control  variables  which  must  be  specified 
on  the  input  deck.   They  are  used  to  switch  various  functions 
in  the  program,  and  to  terminate  the  computations. 


KSTOP   ±1 


0 
JSTOP   +1,  0 

-1 

LSTOP   +1 

0 

-1 


(+1,  no  PICTR) 
Read  Type  1  data  cards (-1,  PICTR  at  ut  =  utl, 

Kt  =  Ktl  for  RAO  for  MAIN  1) 
(-1,  PICTR  of  short  crested 
seas  for  MAIN  2) 

Stop  execution  of  the  program! ! 

Main  1  Program  executed 

Main  2  Program  executed 

Read  Type  2  data  cards 

Stop  execution  of  program! ! 

Read  Type  1  data  cards;  switch  back  to  Main  1 
Program 


Examples  of  sample  data  decks  will  resolve  any  confusion 
caused  by  the  above  statements.   The  intent  of  having  the 
switching  variables  was  so  as  to  be  able  to  run  Main  1,  Main  2, 
or  both  together,  and  at  the  same  time  use  some  of  the 
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subscripted  computer  variables  (matrix  elements)  as  dummy 
holders  for  intermediate  calculation  steps.   Unless  the  data 
decks  examples  are  followed,  there  are  probably  several  non- 
workable  combinations  of  the  switching  variables. 

Data  Card  Input  Deck  Format 

The  general  format  for  all  the  variables  is  floating  point 
decimal.   the  exceptions  are  JSTOP,  KSTOP,  and  LSTOP,  which 
are  fixed  point  decimal  (12) . 

There  are  three  different  types  of  data  cards.   To  run 
Main  1,  types  0,  1,  &  3  are  used.   To  run  Main  2,  types  0,  1, 
2,  &  3  are  needed.   The  data  values  on  types  1  and  2  need  not 
be  identical;  the  reason  for  having  type  2  data  cards  was  to 
permit  a  number  of  variations  of  data  without  retyping  a  whole 
type  1  deck. 
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Card  #  Column 

1-10   11-20   21  -  30   31  -  40   41  -  50   51  -  60   61  -  70   71  -  80 
(column  1-2  KSTOP) 
(column  3-4  JSTOP) 


1A 

IB 

1C 

ID 
IE 


BG 


K 


]iSt 


Dl 


GM 


Y./m 

v   ' 


utl 


D2 


Cb 
At 

D3 


A 


blank 

0)1 

Fl 


B 
Le 

Aoj 
F2 


T 


K./I 

p   x 

F3 


\l 

B 
blank 


II 


1/3 


Ak 

t 

blank 
blank 


/*3 X-axis  label,  colums  1-40 — **f*2 Y-axis  label,  columns  41-80 — &/ 


2A 
2B 
2C 


1/3 
U     H  '  jit       K       (column  41-42  LSTOP)  43-80  blank 


-Title  to  be  printed  below  short  crested  sea  results,  col  1-80 — fe/ 
-X-axis  lable  colums  1-40 — &y^ Y-axis  label,  columns  41-80 £>/ 


-blank  card- 


Combinations  of  data  cards  for  successful  computer  runs: 

1.  Find  optimum  tank  parameters  (Main  1) 

(1A  -  IF)  ,  (1A  -  IF)  , (1A  -  IF)  ,  3;  set  KSTOP  =  01,  JSTOP  =  00 

2.  Find  response  versus  heading  angle  (Main  2) 

(1A  -  IF)  ,  (2A  -  2C)  ,  (2A  -  2C) (1A  -  IF)  ,  (2A  -  2C) (2A  -  2C)  ,3 

Set  KSTOP  =  01  or  -1,  JSTOP  =  01,  LSTOP  =  01  or  -1. 
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7.    Sample  Output  Data 


SIG  WAVE  HT  =   10.0  HEADING  ANGLE  =  9  0.0  DEG 
SHIP  SPEED  =  43.00  EFFECTIVE  SHIP  LENGTH  =  260.0. 
SHIP  NAT  FREQUENCY  =  0.5142  SHIP-TANK  COUPLING  =  0.2  000 
GM=   3.3  60   SHIP-TANK  DECOUPLING  FREQ=   1.04 
LB  =  0.2167  LY=  8.9732  Gl=   3.3155   G2=  11.5796 
KS  =  0.0453  KY=-2. 33446   MAS=    231884.  IX=    75130336. 


KT 

FI1/3 

1 

-0.400 

3.544 

-1.400 

5.783 

-0.900 

4.870 

-0.400 

3.544 

0.100 

4.966 

-0.150 

4.401 

-0.400 

3.544 

-0.650 

4.089 

-0.525 

3.733 

-0.400 

3.544 

-0.275 

3.706 

MT=    1. 

01    Fll/3= 

3.71   ETA1/ 

TANK  LINEAR   DAMPING 

=  -0.27 

KT 

FI1/3 

2 

-0.400 

3.474 

-1.400 

4.202 

-0.900 

3.602 

-0.400 

3.474 

0.100 

5.162 

-0.150 

4.300 

-0.400 

3.474 

-0.650 

3.374 

-0.900 

3.602 

-0.775 

3.472 

-0.650 

3.374 

MT=    1. 

,51    FIl/3= 

3.37  ETA1/3 

TANK  LI 

[NEAR  DAMPING  = 

=  0.65 

5.67 


6.33 


KT       FI1/3  3 

-0.400  7.451 
-1.400  6.939 
-2.400  6,430 
-3.400  6.190 
-4.400  6.067 
-5.400  5.996 
-6.400  5.952 
-7.400      5.922 

MT=    0.51    fIl/3-     5.92  ETAl/3=       0.32 
TANK  LINEAR  DAMPING  =   -7.40 

MI(MIN)+    1.35  FIl/3=  -1.8098 
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LONG  CRESTED  SEAS;   SPEED=   34.00   SIG  WAVE  HT=15.0) 


MT=  KT= 

1.20  -1.00 

WS=  LB= 
0.40      0.22 


LT= 

0.20 
KS= 
-0.06 


HEADING  ANGLE   SIG  ROLL  ANGLE   SIG  TANK  ANGLE  (DEGREES) 


-1 


0.0 

15.0 

30.0 

45.0 

60.0 

75.0 

90.0 

105.0 

120.0 

135.0 

150.0 

165.0 

180.0 


0.00 
2.23 
4.61 
6.92 
4.44 
2.08 
1.29 
0.87 
0.59 
0.39 
0.24 
0.11 
0.00 


0.00 
4.27 
9.34 
15.78 
13.33 
8.28 
5.59 
3.99 
2.83 
1.93 
1.20 
0.58 
0.00 


SHORT  CRESTED  SEAS 


-1 


MT= 

KT= 

LT= 

1.20 

-1.00 

0.20 

WS= 

LB= 

KS= 

0.40 

0.22 

-0.06 

HEADING  ANGLE   SIG  ROLL  ANGLE   SIG  TANK  ANGLE  (DEGREES) 


0.0 

3.41 

7.58 

15.0 

3.42 

7.77 

30.0 

3.44 

8.23 

45.0 

3.40 

8.68 

60.0 

3.19 

8.71 

75.0 

2.75 

8.08 

90.0 

2.14 

6.83 

105.0 

1.48 

5.25 

120.0 

0.93 

3.72 

135.0 

0.57 

2.52 

150.0 

0.37 

1.72 

165.0 

0.26 

1.27 

180.0 

0.23 

1.13 

SHIP  #  2 

HIGH  TANK     LT=0. 

20 
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Appendix  D.   Other  Solution  Methods 

1.  Analog  Solution 

The  equations  of  motion  were  put  into  their  proper  non- 
dimensional  form  for  electrical  analog  simulation.   However, 
it  was  considered  necessary  to  further  investigate  the  numer- 
ical solution  in  order  to  properly  normalize  the  non-linear 
variables,  and  especially  coefficient  expressions.   Using  the 
results  of  this  paper,  an  analog  simulation  is  a  trivial 
problem  if  the  modeling  of  the  coefficients  can  be 
accomplished . 

2 .  Numerical  Integration 

In  order  to  perform  this  technique,  the  three  coupled 
equations  forming  a  fifth-order  system  must  be  algebraically 
broken  down  into  five  coupled  first  order  differential  equa- 
tions.  This  was  done  using  the  set  of  dimensionless  equa- 
tions in  preference  to  the  dimensional  set.   The  five  state 
variables  chosen  to  represent  the  fifth-order  system  were: 
<$> ,    <t>  i      v*,  ri,  and  n .   One  of  the  excitation  terms,  the  iX 
term,  was  dropped  to  make  the  algebra  a  little  more  tractable 
A  parametric  study  was  done  in  order  to  justify  this  step. 
The  below  figure  shows  that  for  the  test  case  investigated 
for  a  "normal"  sized  ship,  the  iX   term  has  practically  no 
bearing  on  the  excitation  f or w  ^>  0.2,  which  was  previously 
shown  to  be  the  lower  boundary  of  the  excitation  frequency 
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range  of  interest.   Therefore,  the  sway  excitation  term  was 
shortened  to  X  \b    and  roll  excitation  to  (1-X,  A  )i|i.   This  kept 
the  excitation  term  uni-dimensional . 


90°  - 


X      =    10. 
for : /   y  > 

o 


1.5 


1.0  Yv      +    Yv 

l  to  CO 


0.5 


CO 


Yv 


CO 


o 


Figure  Dl.   Sway  Excitations. 

The  five  coupled  ordinary  differential  equations,  in 
matrix  form,  were  found  to  be, 


d/dT 


* 

= 

$ 

= 

V* 

• 

0 

= 

n  _ 

= 

VBo      h/^o      V3o      V3o      V3o 


z-j/Yj      z2/Y2      z3/Y2      z4/y2       z5/y2 
6^2      62/a2      63/a2      64/a2      65/a2 


+ 


VPo 


Z6/Y2 
66/a2 


*(t) 


V*o 


z/62 

6/a2 


L    o 


p*  + 


V6o 
0 

VY2 

Va2 
0 


-1/*R2»o" 

0 

Z9/Y2 
69/a2 

• 

PI  Pi 

*3 

0 
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The  above  matrix  equation  is  separated  into  a  ship-tank 
linear  matrix,  and  excitation  matrix,  a  controller  matrix, 
and  a  ship-tank  non-linear  matrix.   The  coefficients  in  the 
matrices  are  all  constants  for  a  given  set  of  conditions.   A 
separate  equation  for  p*  should  be  added  for  a  particular 
contoller . 

The  identities  for  the  29  new  coefficients,  3  ~37,  z,-z  , 


6,-<5q,  and  a  -a_  are  given  below. 


2 

rWi     At    .                A      \ .               ai 

al    =    [ 2 2]               a2    =    [— 2    "    7~T-  ]               a3    =   ^7 

yyt   Y2         yst                                 yt           V   ^2  2 


.     2      2 
Ab   Yl 


3Q   =    [1-   c^   +    Xb   Yl y~   1 


2  "yt   Y2 


,  (X  _  A  — Ot_A     K 

63   "   Ky[Ab(1-^»-   -M      !  »4    =  -4^  65   ■   -*t(1+03' 

'2         y         Yn  t 

yt    '  2 


P6   "    I1+Vy   +  ^T  +    Xbh  87    ■   °3Xt 

"yt   Y2  2 
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34         XtKt 

6i     =     [aiiT  "  "TT"1 


35 

65     =     [aiF"  "  xt] 

o 


6C      =      ^[1    +    A,  A      +    A     (-±-| —  +   -^-±)]    +   -M£ 
6  3  t»yy2  yn  2 


'yt    '2 


yt    '2 


37 
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"alkQ 


8  3 
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69      :   ~~2 


*R    3o 
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13                   b   y  y    ,     a~                     2                 3                        2 
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*yt      2 
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The  matrix  equations  may  be  drawn  as  a  signal  flow  graph 
in  block  form.   The  signal  flow  graph  is  useful  for  inter- 
preting exactly  what  the  basic  relationships  are  between  all 
the  variables.   This  information  cannot  be  extracted  from  the 
simpler  form  of  the  three  coupled  equations. 
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m  "  1 1 1  r  mi  wnr'^TtTwnr, 


V 


■j      v  ,'  w  ,r 


Figure  D2.   Ship  Anti-Rolling  Tank  System 
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Simulation  of  roll  by  the  use  of  the  above  equations 
without  the  necessity  of  directly  using  the  spectral  density 
function  technique  could  be  accomplished  by  taking  an  excita- 
tion made  of  randomly  phased  harmonic  waves,  in  their  proper 
proportions,  and  integrating  for  a  ship  time  duration  of 
enough  cycles,  say  a  thousand,  to  get  a  meaningful  sample  for 

determining  i>         .   In  the  interests  of  amount  of  effort 
rms 

required,  this  should  be  done  subsequently  to  finding  the 
approximate  solution  by  the  complex  algebraic  method  adopted 
by  the  author. 
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